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ABSTRACT 

Using archival SDSS multi-epoch imaging data (Stripe 82), we have searched for flares due to the 
tidal disruption of stars by super-massive black holes in non-active galaxies. Two candidate tidal 
disruption events (TDEs) are identified, using a pipeline with high rejection efficiency for supcrnovac 
(SNe) and Active Galactic Nucleus (AGN) flares, and minimal selection bias. The flares have optical 
black-body temperatures 2 x 10 4 K and their cooling rates are very low; their observed peak lumi- 
nosities are M g = —18.3 and —20.4 {yL v = 5 x 10 42 , 4 x 10 43 ergs -1 , in the rest-frame), qualitatively 
consistent with expectations for tidal disruption flares. The properties of the TDE candidates are 
examined using i) SDSS imaging to compare them to other flares observed in the search, ii) UV emis- 
sion measured by GALEX and Hi) spectra of the hosts and of one of the flares. Our pipeline excludes 
optically identifiable AGN hosts, and our variability monitoring over 9 years provides strong evidence 
that these are not flares in hidden AGNs: 1) the luminosity increases in the candidate TDE flares are 
estimated to be greater by factors 4 and 15 than in any of the variable AGNs monitored, and 2) the 
TDE candidates' hosts are much quieter in the seasons not including the primary flare than are hosts 
of AGN flares (less activity on average in each of the non-primary-flare seasons than in 95% of AGNs, 
with a cumulative probability estimated to be < 1CU 5 for a flaring AGN to have all other seasons as 
quiet as for the TDE candidates). The spectra and color evolution of the flares are unlike any SN 
observed to date, and their strong late-time UV emission is particularly distinctive. These features, 
along with the high resolution with which they are placed at the nucleus, argue against these being 
first cases of a previously-unobserved class of SNe or more extreme examples of known SN types. 
Taken together, the observed properties are difficult to reconcile with a SN or AGN-flare explanation, 
although an entirely new process specific to the inner few-hundred parsecs of non-active galaxies can- 
not be excluded. Based on our observed rate, we infer that hundreds or thousands of TDEs will be 
present in current and next-generation optical synoptic surveys. Using the approach outlined here, 
a TDE candidate sample with 0(1) purity can be selected using geometric resolution and host and 
flare color alone, demonstrating that a campaign to create a large sample of tidal disruption events, 
with immediate and detailed multi-wavelength follow-up, is feasible. A by-product of this work is 
quantification of the power-spectrum of extreme flares in AGNs. 



1. INTRODUCTION 

When a star passes too close to a super-massive black 
hole, the tidal shear overcomes the star's self-gravity and 
the star is disrupted. Much of the stellar debris is ejected 
from the system, but some fraction remains bound to the 
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black hole and is accreted , resulting in a week- to year- 
long electromagnetic flare (jRecs 1988). The fall-bac k rate 
is exp ected to follow a power- law of index —5/3 (iRees 
Il98ct iPhinnevI 119891) . but see lLodato. King, fe Pringle 
(|2009|) for predicted deviations from this canonical scal- 
ing; the light -curve does not in general show the same be- 
havio r, (e.g.. iStrubbe fc Quataertll2009t lLodato fe Rossil 
120 111) . For M BH < 10 7 M Q the initial fall-back rate is 
super-Eddington and the emission is usually assumed 
to have a black body spectrum. The predicted tem- 
peraturcs range from ~ 10 4 K (jLoeb fc Ulmerl I1997D 
to, more commonly, ~ 10 5 K (e.g., lUlmerl Il9990 . A 
radiativcly driven wind, existing as a consequence of 
the super-Eddington fallback, may dominate the emis- 
sion of the tidal disruption event (TDE) at early times 
(|Strubbe fc Quataertl|2009l ). Model predictions generally 
depend on parameters which are quite uncertain, so ob- 
servational input is needed to constrain the modeling. 

Detections of tidal disruption events are of interest 
for a number of reasons, including: 1) The light emit- 
ted after the disruption depends sensitively on the black 
hole mass and spin, hence a large sample of TDEs 
will allow properties of back holes to be studied with- 
out relying on scaling relations with global parameters 



of ga l axies (iGebhardt et all 120001: [Ferrarese fc MerritH 
[20001 IGraham et al.1 120011 IMarconi fc Hunt! 120031) . 2) 
TDEs may be our only probe to obtain a large sam- 
ple of dormant super-massive black holes (jFrank fc Reed 
fl97llLidskii~fcO zernoillT979f ). 3) A particularly intrigu- 
ing application is testing the existence of intermediate 
mas s black holes in globular cluste rs and dwarf galax- 
ies (jRamirez- Ruiz fc Rosswol[2009l ). 4) For black holes 
with mass Mbh ^ 1O 8 M , the tida l disruption radius 
lies inside the Schwarzschild radius ([Hillsl Il975l ) . hence 
a TDE survey that covers a sample of galaxies with a 
wide enough central black hole mass range is in prin- 
ciple sensitive to whether super-massive black holes do 
have an event horizon. 5) Detailed observations of the 
emission from a large sample of tidal disruption events 
will provide a new arena for testing our understanding 
of accretion physics and may constrain properties of the 
disrupted stars. 

A number o f candidate TDEs have been identified in 
X-ray surveys (IBade et al.lll996l)Komossa fc Badelll999l: 
Donley et all 120021: lEsauei et al.l 120081 ICappelluti et al.l 



20091: iMaksvm et al.ll2010l ) — for a review see iKomoss 



2001— and in the UV GALEX Deep Imaging Survey 



Gezari et al.l 120061 [20081 [200 9b). However establishing 



that a candidate tidal flare found in UV and X-ray sur- 
veys is not an exceptionally variable AGN is hampered 
because, although the amplitude of AGN variability at 
these wavelengths is much greater t han in the optical 
(|Maoz et al.l 120051; iSaxtonet ail 120111 ). the range of vari- 
ation has not yet been well-characterized. With the ad- 
vent of optical transien t surveys, such as the Palomar 
Transient Factory (PTF: lLaw et al.ll2009[ ). Pan-STARRS 
(|Chambersll2007| ) . Catalina Real-Time Transient Survey 
(CRTS) and later the La rge Synoptic Survey Telescope 
(LSST: uvezlc et al.l [2001 . there should be many candi- 
date TDEs. The challenge is to eliminate the far more 
common flares from variable AGNs and SNe to a) con- 
vincingly exclude these backgrounds as explanations for 
individual events and b) efficiently produce a high-purity 
sample of TDE candidates such that expending resources 
following up uninteresting events can be reduced to an 
acceptable level. 

Until now, the feasibility of identifying probable tidal 
flares with an optical transient survey alone has not been 
demonstrated. Here we present two candidate stellar 
tidal disruption events detected in archival multi-epoch 
imaging data of the Sloan Digital Sky Survey (SDSS), 
showing that the difficulties of identifying such events in 
an optical survey can be resolved. There are several im- 
portant advantages of using the SDSS data. First is the 
very large sample of galaxies, with many having spec- 
tra; this sample enables one to classify flares into well- 
defined categories, such that the validity of our method- 
ology can be demonstrated. Second, the spatial resolu- 
tion of SDSS is adequate to exclude a very high fraction 
of SNe from the nuclear-flare sample, allowing SNe to be 
rejected without imposing cuts based on flare properties. 
Third, SDSS observed the galaxies in Stripe 82 over typi- 
cally 7 seasons with a mean of 70 observations in all; this 
allows hidden AGNs to be rejected based on variability 
in the non-peak seasons. Having observations in three 
or more filters is also useful, providing color information 
that is valuable in confirming that the flares are not due 



to AGNs or known types of SNe. 

The selection pipeline we employ reduces the back- 
ground from SNe and variable AGNs by two orders of 
magnitude or more, but positive determination that the 
two flares which pass the pipeline are in fact TDEs re- 
quires detailed comparison between the properties of the 
flares and those of AGNs and SNe. This we do with 
multiple tools. We have obtained spectra of both hosts 
and of one of the flares, found archival GALEX post- 
flare observations of both host galaxies and a pre-flare 
observation of one of them, and archival Catalina Real- 
Time Transient Survey data to extend the light curve of 
one of the flares to earlier and later times than observed 
by SDSS. Analysis of these observations and comparison 
between the properties of the TDE candidates and SNe 
and AGNs is reported. 

In addition to demonstrating the feasibility of 
optically-initiated TDE surveys, discovery of these two 
flares gives needed observational insight into the tidal 
disruption phenomenon. An important consequence of 
our method of detection - which does not rely on flare 
properties beyond requiring them to be nuclear to reject 
SNe and uses host properties alone to reject AGNs - is 
that selection bias is minimized. With just two events 
we have only begun to scratch the surface, but the prop- 
erties of the flares can already test and inform theory. 
Our candidate TDEs are b oth much more luminous i n 
the optical than predicted (jStrubbe fc Quataertl |2009| ) . 
although with adjusted parameter choices the fit can be 
improved and some ingredients may still be missing in 
these early, simple models. 

The outline of the remainder of this paper is as fol- 
lows. In Section [2] we describe our pipeline for TDE se- 
lection. In Section [3] we present the observed properties 
of the final products of the pipeline: two candidate tidal 
disruption flares. Follow-up observations with other in- 
struments are reported. A detailed comparison to flares 
of active galaxies and supernovae is given in Section [U 
leading to the conclusion that the two TDE candidates 
are indeed likely to be tidal disruption events. Section 
[5] compares our TDEs to candidates reported in other 
wavelengths and compares their observed properties to 
theory. The implications of this work for detecting TDEs 
and obtaining a relatively pure sample of them in future 
optical surveys is discussed in Section [6] We close with 
a summary (Section [7]). 

All mag nitudes are quoted in the AB system 
Ol^lT97l an d are corrected for Galactic extinction 



Schlegel et al.l [l998h . We adopt a standard cosmology 
with H = 72kms- 1 Mpc" 1 , n m = 0.3 and n A = 0.7. 



2. IDENTIFICATION OF TIDAL FLARES 

The flowchart in Fig. [1] summarizes our search graph- 
ically. In the following sections we discuss the steps in 
this chart in more detail; for additional details see S. 
van Velzen, G. R. Farrar et al. (2011). To achieve the 
goal of obtaining a large and uniformly selected sample 
of tidal disruption events (TDEs), it is necessary to ac- 
curately identify their flares in large volumes of optical 
data. Here we use SDSS observations of 2.5 x 10 6 galax- 
ies in "Stripe 82" , that have been observed typically 70 
times over a nine year baseline. In Section ^. ll wc explain 
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Summary of catalog cuts. 
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Hidden 
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Fig. 1. — Flow chart summarizing the flare and TDE identifica- 
tion pipeline, discussed in Section[2] The numbers in the grey boxes 
indicate the number of objects in ea ch class. The first three steps 
are discussed in Sections l2,l.Hf2,2.2l the last three step s summarize 
the selection for TDE candidates discussed in Sections 12. 2. 3lf2. 3. 31 

how a sample of 342 flares was extracted from this data 
set. Supernovae or other stellar flares are not of interest 
in a search for TDEs, and in Section |2~2"1 we explain how 
they are removed by discarding flares offset from the cen- 
ter of the host. The final sample for our study consists 
of 42 nuclear flares having more than two observations 
in the flaring state. 

The rejection of active galactic nuclei (AGN) is dis- 
cussed in Section 12.31 Eliminating host galaxies with 
spectroscopic and photometrically identified AGNs re- 
moves all but 5 nuclear flares. The critical challenge is 
to identify and exclude host galaxies with active nuclei 
which arc too weak or heavily obscured to be identified 
by usual spectroscopic and photometric criteria. We can 
identify variable unrecognized AGNs directly by their 
variability, thanks to SDSS's multi-year monitoring of 
the hosts. Three of the five candidates do display vari- 
ability in other seasons, consistent with the variability 
we measure in our identified AGNs, while two others do 
not, at the < 10~ 5 CL level (see Section |4~T|) . On ge- 
ometric grounds alone, the probability these two flares 
are SNe is < 0.5% (see Section |4~2|) . assuming the distri- 
bution of SNe follows the stellar light. Thus these two 
flares are strong candidates to be stellar tidal disruption 
events and we refer to them as TDE1 and TDE2 below. 

2.1. Selecting flares in SDSS observations 

The ~ 300 deg 2 Sloan Digital Sky Survey (SDSS; 
lYork et all 120001) multi-epoch imaging data Stripe 82 
dSesar et all 120071: iBramich et al.1 120081 : iFrieman etafl 
l2008HAbazaiian et al\\2(M ) was the starting point of our 
search. Using the SPSS morphological star-galaxy sep- 
aration (jLupton et alj|2001t iStoughton et all 120021) and 
the s tandard checks on quality flags (|Stoughton et al.l 
we extracted ~ 2.5 10 6 galaxies from this data 
set. The typical number of observations per galaxy is 70. 

2.1.1. Catalog cuts 

In order to focus our analysis resources and time most 
efficiently, we first selected candidate flares from the cat- 
aloged parameters, which reduced the number of galaxies 
by two orders of magnitude. Galaxies were required to 
have m < 22.5 in at least 3 of the 5 SPSS bands (u,g,r,i,z; 
iFukugita et "all H99l iSmith et al.l [2002). A minimum 
flux increase of 10%, measured at the 7-a level, is the 
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identified as extended in co-add 
X 2 /DOF > 5 



- peak 



> 1.1 



(Fpeak ^baseline) /<Tp ca k > 7 
(-^peak -^baseline)/ fT rms > 3 



> 3 
all 

g,r, or i 

> 2 

> 2 
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Note. — The first two cuts are designed to 
obtain a clean, flux limited sample of galax- 
ies. The co-add runs are obtained by ad- 
dition of nearly all S tripe 82 imaging data 
IjAbazaiian et al.ll2009h . J. Annis et al. (2009), 
in preparation. The goal of the remaining cuts 
is to select flares. For all cuts we use the Pet- 
rosian flux UStoughton et al.ll2002T) . The sub- 
script 'mean' is the inverse-variance-weightcd 
mean using all observations, while the subscript 
'baseline' refers to the non-flare observations 
only, x 2 i s calculated using mean flux as a 
model for the galaxy light curve. These cuts se- 
lect ~ 2x 10 4 candidate flares from the ~ 2x 10 6 
galaxies with m < 22.5 in Stripe 82. 

most important requirement that was imposed to find 
flares. The cuts on the cataloged parameters were chosen 
to be "soft" (i.e., high efficiency, low purity) and yielded 
21,383 potential flares. See Table[T]for a summary of the 
catalog cuts. 

2.1.2. Difference imaging 

Next, with this data set that is a factor of 100 smaller 
than the starting one, we applied a more rigorous anal- 
ysis method: image subtraction. Given that there are 
typically 70 observations of a galaxy, we have many im- 
ages of the host of the flare; from these we selected the 
observations with the best seeing and lowest sky level 
to be the reference images for subtraction. First we 
implemented a simple and relatively fast "direct sub- 
traction" method, which allows quick determination that 
8834 flare candidates are spurious because they show no 
flux in the difference image. For each of the remaining 
galaxies hosting a flare, nine reference images (of size 
1' x 1') were cro ss convolved using a mo dified version of 
the software by lYuan &; AkerloJ (|2008l ) and subtracted 
to obtain nine difference images, for each filter in each 
night. The convolution kernels were determined using 
nine reference stars that were selected to lie close to the 
host galaxy. The mean flux of the nine difference images 
and its standard deviation were computed for each pixel, 
after applying a clipping algorithm to reject pixels that 
lie more than 3<7 from the median. The flux, F, in the 
mean difference image was computed using an aperture 
of 2 times the full-width-half-maximum (FWHM) of the 
point spread function (PSF) measured in the mean im- 
age of the convolved reference stars. The minimum flux 
in the difference image was required to be m < 22. In 
addition, we require for each band that the flare is de- 
tected at the 7-a level. We further require that the shape 
of the source in each difference image is well-fit by the 
PSF extracted from the mean convolved reference stars 
for that image. 

Of the flares meeting the above conditions, 583 satisfy 
the requirement of detections in at least two bands. Good 
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agreement on the flux in the diffe rence image is obtained 
for the light curves produced by iHoltzman et all (12008D 
of su p ernovae from the SPSS SN Survey dFrieman et al.l 
120081 iSako et al.1 120081 ) which were detected indepen- 
dently by our pipeline. 

2.f .3. Rejecting moving objects 

Solar system objects, when seen in front of a galaxy, 
can fake a flare. SDSS checks for moving objects by 
measu ring the position of the centroid across different 
filters (jlvezic et al.ll2001[). We eliminate moving objects 
flagged by SDSS (|Stoughton et al.ll2002t ). and identify 32 
additional solar system objects using the position of the 
centroid in the difference image. By requiring detections 
in the difference image in at least two observing nights, 
false flares due to solar system objects are eliminated and 
we are left with 41 9 flares. 

2.1.4. Manual rejection 

We remove galaxies with a spectroscopic redshift > 1.2, 
because they are mis-classified as extended. This re- 
quirement removes 7 quasars from the sample. Finally, 
we inspect the difference images of the remaining 412 
flares to search for bad subtractions or other anoma- 
lies. The number of flares rejected after this inspection 
is four, which is quite low considering the large number 
of galaxies processed by the image subtraction pipeline 
(~ 2 x 10 4 ). In Fig. [5] we show examples of rejected 
and good subtractions. Additional quality cuts on the 
distance between the flare and the center of its host are 
discussed below. 

2.2. Selection of nuclear flares 

An accurate measurement of the distance between the 
center of the host galaxy and the flare is crucial for ob- 
taining a clean nuclear flare sample (i.e., removing SNe). 
In Section l2.2.1l we discuss the details for accurately mea- 
suring this distance and in Section 12.2.21 we explain how 
the host-flare distance is used to define the nuclear flare 
sample. 

2.2.1. Precise determination of host-flare distance 

In this section wc consider first the accuracy with 
which we determine the position of the flare with respect 
to the Sersic center of the host, and then the extent to 
which the Sersic center of the host is an accurate deter- 
mination of the true galactic center. 

The location of the flare relative to the host is found 
by fitting a point source corrected by the PSF to the 
difference image. We repeat these measurements for all 
reference images in all nights with detections in the g, 
r or i bands; this yields at least 9 x 2 x 2 = 36 mea- 
surements of the distance d between the flare and the 
center of the host. The uncertainty on this distance in 
pixel coordinates, o~ x and o~ y , is obtained from the stan- 
dard deviation of these measurements. As seen in Fig. 
[31 the typical measurement uncertainty on the distance 
between the host and the flare is 0.1 pixel. We use this 
as a lower limit for cr x ,y To divide flares into categories 
of clearly nuclear, clearly non-nuclear or ambiguous, we 
define cr^ = ^J(xo x ) 2 + {ya y ) 2 /d, which is a measure of 
the uncertainty on the host-flare distance. 




(a) Rejected: bad subtraction, one noisy filter 









(b) Rejected: bad subtraction, peak displaced between niters 




(c) Good subtraction: near the flux limit 




(d) Good subtraction: off-center flare (SN) 

Fig. 2. — Examples of bad and good image subtractions. Image 
size is 1' X 1'. Left to right: flare image, mean reference image 
and difference image. The bad subtractions, images (a) and (b), 
were rejected by manual inspection; image (c) shows the quality of 
a subtraction that is relatively faint, m r = 21.5 in the difference 
image. Image (d) shows an example of an off-center flare detected 
at r = 0'.'4 ± 0'.'02. 



In the above, d is defined with respect to the coordi- 
nates of the optical center of the host obtained by fitting 
a Sersic profile to the convolved reference image. How- 
ever a deviation from azimuthal symmetry in the galaxy 
(e.g., an H II region) may cause the center found by fit- 
ting a Sersic profile to be offset from the true center, thus 
causing an error in determining the separation of the flare 
from the galaxy center. We have investigated the impact 
of this on our accuracy, galaxy-by-galaxy, by comparing 
the center obtained from a Sersic fit to the center ob- 
tain by fitting the PSF to the galaxy. To enhance any 
deviation, we subtract 50% of the best-fit azimuthally 
symmetric model galaxy from the original image before 
fitting the PSF. As shown in Fig. |4j the two methods 
on average yield the same answer. Visual inspection of 
the few galaxies with large (> 0.1 pixel, i.e., 0'.'04) de- 
viations between the two ways of measuring the center 
shows that they indeed have more complicated shapes 
or consist of close pairs of galaxies. In those rare cases 



Optical discovery of tidal disruption flares 



5 



40 t 



30 r 



20 r 



10 r 



x [mean, a]: 0.021, 0.122 
y [mean, <rj: -0.022, 0.095 



-0.4 -0.2 0.0 0.2 0.4 
distance between host and flare (pixel) 

Fig. 3. — Histogram of the distance between the center of the 
host and the flare for all 186 nuclear flares (eq. [T}. The dashed 
lines display the best-fit Gaussian distribution. As expected for 
nuclear flares, the distribution peaks at zero. The mean accuracy 
on the distance between the host and the flare for the nuclear 
sample is about 0.15 pixel or 0'.'06. This is similar to the SDSS 
single ep och astrometric a ccuracy for point sources brighter than 
m r ~ 20 <Pier et alj|2003ft . 
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Fig. 4. — The distance between the center of the galaxy obtained 
by fitting a Sersic profile and the center obtained by fitting a PSF. 
Irregularities in the galaxy may push the center obtained by the 
Sersic fit away from the true geometrical center. To enhance this 
potential systematic uncertainty, we subtract 50% of the best-fit 
Sersic model galaxy from the original image before fitting the PSF 
to find the center. We find that for almost all galaxies the two 
methods of locating the center agree very well: the Sersic profile is 
not significantly affected by irregularities. (1 pixel corresponds to 
0'.'4) 



that the magnitude of the difference between the PSF 
and Sersic centers in either coordinate is larger than the 
uncertainty on the separation between the flare and Ser- 
sic center, we take the flare positional uncertainty cr x ,y, 
to be the former. 

Now, with the positional uncertainties understood, we 
proceed to define our flare sample. First, we set the 
maximum distance between the host and the flare to d < 
1" and remove all transients detected beyond this radius 
from our flare sample. In addition, we demand that the 
distance between the host and the flare is measured with 
a minimum accuracy of ad < O'.'l. After applying these 
cuts, we are left with flares in 342 galaxies. 

2.2.2. Separating nuclear flares from non-nuclear flares 
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Stellar PDF 
Total PDF 








0.0 



0.2 



0.4 0.6 
d (arcsec) 



0.8 



1.0 



Fig. 5. — The observed distribution of host-flarc distances (thin 
black histogram, with error bars given by Poisson statistics), the fit 
(filled black histogram) with the decomposition into nuclear flares 
(blue histogram) and SNe (purple histogram), all normalized to 
unit area. The only free parameter in the fit is the fraction of 
stellar-distributed flares, P(SN). The top panel shows the d dis- 
tribution for the full sample of 342 flares; the bottom panel shows 
the flare sample that remains after requiring at least 2 observations 
after the peak of the flare, imposed for TDE analysis. 



In this section we use the distance d between the center 
of the host galaxy and the flare to divide the sample of 
342 flares into two well-separated samples of 186 nuclear 
flares and 85 off-center flares, plus 71 flares that are not 
used because their classification is ambiguous. 

We begin by determining the overall fraction of SNc 
in the full sample of 342 flares, P(SN). To do this wc 
model the distribution in d as a sum of nuclear flares 
and stellar-distributed flares (i.e., SNe). The unsmeared 
probability density function (PDF) for nuclear events is 
a delta function at d = 0. We assume that stellar flares 
trace the stellar light, as is justified by the discussion in 
Section l4~2l The flare detection efficiency as a function 
of d has been modeled and tested against observations, 
as will be reported in detail elsewhere; for the purposes 
of this paper it is sufficient to note that the detection effi- 
ciency should be at most weakly dependent on d, except 
for a possible reduction in efficiency near the nucleus. 
We confirm this a posteriori below. 

Under these assumptions, the unsmeared PDF for 
stellar-distributed flares at distance d + Ad from the host 
center is given by the sum total galaxy flux in this inter- 
val. For each host galaxy, the surface brightness profile, 
F(d), is taken to be that of the PSF-corrected model 
galaxy corresponding to the best-fit Sersic parameters in 
the r-band. We smear the nuclear and stellar PDFs of 
each galaxy according to the measured uncertainty in the 
flare separation for that galaxy, a x and o~ y . 

We make an unbinned maximum likelihood fit to 
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the sum of the resulting stellar-distributed PDF, times 
P(SN), plus the nuclear PDF, times l-P(SN), to the ob- 
served d distribution, with the result P(SN) = 0.34±0.04 
(90% CL). The observed and predicted distributions are 
shown in Fig. [5] (top panel). The quality of this fit, in 
which the SN fraction P(SN) is the only free parameter, 
validates the assumption that the detection efficiency is 
independent of d for d > 0'.'2 and shows that we under- 
stand the distribution of host-flare separations. 

Based on Fig. [5j we formulate the cut for nuclear 
flares: 

cut = (d/a d < 2) and (d < 0'.'2) . (1) 

This gives 186 on-center flares, with a loss of 15% of 
the true nuclear flares. The probability of a background 
event (i.e., a stellar-distributed flare within this cut), 
can be computed directly from the e?-histogram and is 
P(SN|cut) = 0.052, giving an expected number of back- 
ground events in the total nuclear flare sample of 10. 
This result is independent of whether the <?, r or i-filter 
is used to obtain the surface brightness profile of the 
galaxies. This procedure gives an upper bound on the 
contamination of the SNe in the nuclear sample, since 
a possible reduction in detection efficiency for d < 0'.'2 
would result in a lower predicted number of SNe in the 
nuclear sample (and larger rate of nuclear flares relative 
to SNe, generating the observed flare sample). 

A high-purity sample of flares which are clearly off- 
center from their host is obtained by requiring 

(d/a d > 3) and (d > 0'.'2) . (2) 

This selects 85 flares, of which 27% have been spectro- 
scopically identified as SNe by the Sloan Digita l Sky 
Survey-II Supernova Survey ([Frieman et al J 12008) . We 
assume that the remaining off-center flares are SNe as 
well and take these 85 flares as our SN sample below. It is 
important to note that - because this sample is obtained 
by a cut on the host-flare distance only - the properties 
of the SNe in this sample are not subject to the selec- 
tion biases that exist for spectroscopic SNe surveys, and 
are representative of the properties of SNe that appear 
in the nuclear sample (to the extent that SNe properties 
are independent of their location in the host galaxy). 

2.2.3. Potential TDE sample 

To be able to obtain a well-measured decay rate - a 
valuable diagnostic in separating SNe, AGN flares, and 
TDE flares - we require at least two detections after the 
peak of the flare in the u, g and r bands. (At this step, 
the nuclear - potential TDE - sample is reduced to 42 
flares and SNe sample to 12 flares.) Repeating the d- 
distribution analysis of Section [2.2.2[ the expected num- 
ber of residual SNe in the nuclear flare sample is 0.9. 
The corresponding histogram is shown in Fig. [5] (bottom 
panel). 

2.3. AGN rejection 

AGNs are well-known to be variable and we expect 
that the majority of the 42 nuclear flares that remain 
at this stage of the analysis originate from active black 
holes. Fortunately, AGN hosts are readily identified us- 
ing SDSS spectra of the host galaxy and, with a small loss 
in efficiency, can be identified photometrically too. Af- 
ter applying these methods to identify AGNs in sections 



1.5 : 




Fig. 6. — Color-color diagram of the PSF flux (Stouehton ct al. 
[200% for the 2.5 X 10 6 extended objects in Stripe 82 that were 
used in this work (th in blue line) and the SDSS QSO sample 
HSchneider et all 120071 ) with z < 1. For both samples, contours 
encompass 50, 80 and 90%. The TDE candidates that fall outside 
the QSO locus but are rejected based on additional variability arc 
labelled R1-R3; TDE1 and TDE2 are labelled Tl and T2, respec- 
tively. 



l2~3~T1 and l2~3~2l five potential TDE remain. Motivated by 
previous work on Stripe 82 showing that at least 90% of 
unresolved, spectroscopically confi rmed quasars are y ari- 
able at the 0.03 mag level (rms) (jSesar eFaL|[2003), 111 
section 12.3.31 we use variability beyond the flare season 
to reject three more AGNs, leaving only two flares, we 
label these TDE1 and TDE2. 

2.3.1. Spectroscopic AGN Identification 

SDSS spectra with a median S/N greater than 3 are 
available for about 2/3 of the hosts of the 186 nuclear 
flares. The Princeton ReductiontQ classify galaxies with 
spectra into three classes using a Principal Component 
Analysis: STAR, GALAXY or QSO. Here, QSO does not 
refer to the classical Mb < — 21 luminosity cut, but im- 
plies that the spectral energy distribution (SED) is domi- 
nated by an AGN-like spectrum inside the fiber (3") that 
was used to obtained the spectrum. Of the 42 nuclear 
flares suitable for our TDE analysis, 32 are in hosts that 
are classified as QSOs; we eliminate these from our candi- 
date sample. We also reject hosts in the class GALAXY 
which we can identify as Seyfert galaxies based on either 
of the following criteria: 

• Galaxies that show broad (> 200 km/s line-width, 
at the 7-a level) Ha, H/3, [OIII]A5007, [OII]A3727, 
or MgIIA2799 emission lines. We only consider 
lines detected at the 3-er level with rest-frame 
equivalent width > 5 A. 

• Galaxies that can be classified as Seyf erts using the 
Baldwin. Phillips, fc Terle^ichl (Il98ll ) dia gram: 

/ [OIII]A5007 \ 0M 

H/3 J > log(NIIA6583/Ha) - 0.05 + 

(3) 

()Kauffmann et all I2003D . We apply this formula 
only if all four emission lines are measured at the 
3- a level. 

12 http://spectro.princeton.edu 
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Fi g. 7. — Relative photometry light curves, using the Pctrosian g-band flux, of the hosts of potential TDEs and one Seyfert (see Section 
I4.lt . Each observation is scaled using the flux of nearby (distance < 20') galaxies that are of similar magnitude (Am < 1). The lower 
portion of each panel displays the light curve of one of the ten test galaxies (selected to have Am < 0.3). The observations in the season 
that contained the flare (indicated with red boxes) are excluded in the calculation of x 2 .Variability in the non-flare season is evident for 
all galaxies expect the host of TDE1 and TDE2. 



This eliminates two more objects from the TDE analysis, 
reducing the sample of potential TDEs to 8 flares. 



2.3.2. Photometric AGN Identification 
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We use the photometric properties of the hosts of the 
8 flares that remain at this point to test if they are 
QSOs. We define a locus that contains 9 0% of the 20,710 
spect roscopic QSOs identified by SDSS ([Schneider et al.1 
120071 ) with z < 1 and reject all hosts that fall inside; 
23% of the ~ 2.5 x 10 6 galaxies that were in the original 
sample and processed in the flare search fall inside this 
locus. This cut on QSO locus removes 3 flares from the 
potential TDE sample, see Fig. [51 

2.3.3. Additional AGN variability 

Figure [6] shows that the hosts of five flares in the po- 
tential TDE sample fall outside the QSO locus. We in- 
vestigated the light curves of these hosts to look for ad- 
ditional variability beyond the season that contains the 
flare. Having an efficient but accurate means to quantify 
the flux variability is applicable to a number of studies, 
so in this subsection we give details of the method we use 
to produce relative flux light curves. 

As explained in section 12.1.11 we use the SDSS Pet- 
rosian flux to search for flares in a galaxy light curve. 
The Petrosian flux is designed to yield a robust esti- 
mate of the total flux of a galaxy (jBlanton et al.l [2001t 
iStoughton et al.l [20021) . yet night to night differences in 
seeing or other observational conditions introduce a small 
amount of jitter to our light curves. In most cases, nearby 
and similar galaxies suffer the same attenuation (e.g., 
worse seeing causes a lower Petrosian flux for all faint 
galaxies in the field), hence this jitter can be reduced by 
normalizing the observations of each night using a set of 
reference galaxies. For each host of a flare, we select this 
set by looking for nearby (within 20') galaxies, whose 
apparent magnitude difference (in a given filter) with re- 
spect to the target host is less than one, and for which 
the rms variance in the magnitude in all observations is 
smaller than 30%; this typically yields 70-100 galaxies 
for each host. For each point in the light curve of the 
host, we determine an overall flux correction by averag- 
ing the rcscalings required to bring the flux of each ref- 
erence galaxy to equality with its mean flux. The above 
procedure can be applied to any filter; for our variability 
studies we use the Petrosian g-band flux. 

To test the quality of the resulting light curves, we se- 
lect ten test galaxies from among the reference galaxies, 
chosen to be even more similar to the host of the flare 
by requiring |Ato| < 0.3. Sometimes, fewer than ten test 
galaxies are found; if none are found, no relative flux 
light curve is produced for this host. SDSS provides er- 
rors - let us call them ctsdssm ~ for each observation i 
of each test galaxy t. We test empirically whether these 
error bars are a good representation of the flux uncer- 
tainty after rescaling, by computing Xt > the reduced chi- 
squared assuming a constant flux of the light curve of 
the each test galaxy. If all of the uncertainty in the flux 
measurement were captured in the ctsdssmj the Xt va l~ 
ues would obey a x 2 distribution but they do not. We 
find that the errors reported by SDSS sometimes under- 
predict the true uncertainty of the flux. We therefore 
use the fact that the expectation value of x 2 P er degree 
of freedom should equal 1 for the test galaxies, to ob- 
tain a better estimate of the measurement uncertainty 
on the relative flux the host or any of its test galaxies, 

e.g., O l = CTSDSSiV < Xt >■ 
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Fig. 8. — Histogram of x 2 1 DOF obtained for the relative pho- 
tometry light curves fit to a constant flux, showing the difference 
between QSOs and ordinary galaxies. We use the non-flaring sea- 
sons fo r the Q SO, and all seasons for the nearby test galaxies (see 
Section |2. 3. 311 . For three spectroscopically identified Scyfcrt galax- 
ies with flares we find: x 2 /DOF = {5.14, 2.51, 1.54}. The lowest 
value in this list is larger than 99% of the test galaxies. We may 
thus conclude that flares from Seyfcrt galaxies can also be identified 
based on additional variability of their host. 

The relative flux light curves produced by this method 
using the Petrosian flux in the g-band for the five poten- 
tial TDEs are displayed in Fig. [7] Three of the hosts 
clearly show additional flaring activity: x 2 P er degree 
of freedom of the flux calculated from their relative flux 
light curves excluding the season with the primary flare is 
> 10; for reference, we label these R1-R3. However two 
hosts show no additional variability, with each having 
X 2 /DOF = 1.05. We show in Fig. U that this is sig- 
nificantly smaller than measured for any QSO or Seyfcrt 
galaxy in our sample. These are the hosts of the flares we 
designate as TDE1 and TDE2. In Section |4~T1 we use the 
observations of QSOs and Seyferts to quantify the prob- 
ability that TDE1 or TDE2 are flares from AGN which 
had quiet years in the other observing seasons. 



3. PROPERTIES OF THE IDENTIFIED TIDAL 
FLARES 

In this section we report on the properties of the two 
stellar tidal disruption flare candidates, based on Stripe 
82 SDSS imaging data and observations using other in- 
struments. Although we did not rely on properties of 
the flaring state to select TDEs, to avoid unnecessar- 
ily biasing the selection and because theoretical predic- 
tions are uncertain, the TDEs identified by our pipeline 
prove to be quite distinct from SNe and flaring AGNs. 
This increases confidence that TDE1,2 are not examples 
of familiar phenomena occurring in improbable circum- 
stances, such as SNe accidentally close to the nucleus or 
AGNs which flare dramatically in the midst of a multi- 
year quiet phase. 

We begin with an overview of the two events by giv- 
ing their light curves, cooling rate and flare colors based 
mainly on the SDSS Stripe 82 observations. The SDSS 
images of the TDEs and their host galaxies are given in 
Fig. [9] Tables [2] & [3] summarize the properties of the 
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TABLE 2 
Properties of the TDEs. 



Name 


2 


M a 


L 3 

(xlO 43 erg/s) 


dlnLg/dt 
(xlO^day- 1 ) 


T 

(xlO 4 K) 


d\nT/dt 

(xl0" 3 day~ 1 ) 


d 

(arcsec) 


<%8 
(arcsec) 


TDEl 
TDE2 


0.136 ±0.001 
0.251 ±0.002 


-18.3 ±0.04 
-20.4 ±0.05 


0.54 ±0.02 
4.1 ±0.2 


-1.7 ±0.1 
-0.8 ± 0.1 


2 4 +0.3 

i o 9 +0.07 
1,0 -0.06 


-2 ±4 
-3 ±2 


0.058 
0.068 


0.124 
0.075 



Note. — The rest-frame g-band peak observed absolute magnitudes, luminosities and black body temperatures of the flares, 
measured in their SDSS difference images. We exploited the absence of significant color evolution to improve the accuracy 
of the black-body temperature determination by using the mean flux in each band over the SDSS light curve. In the eighth 
column, d denotes the distance between the center of the host and the flare (see Section [2]l and dgg the 68% confidence radius. 
The host properties (including coordinates) are given in Table [3] 




(b) TDE2 -200 200 400 600 800 

MJD - 53974 



Fig. 9.— SDSS images (T X 1') of the TDE flares and their 
host galaxies. Left to right: flare image, mean reference image 
and difference image for TDEl and TDE2. We see that the hosts 
can be classified as E/S0. The difference image of TDEl shows 
a subtraction artifact at the location of the bright point source, 
which is not subtracted perfectly because our difference imaging 
method is optimized for the center of the field. 

TDEs and their hosts, respectively. Then we describe 
the observations of TDEl and TDE2 obtained with other 
telescopes. Finally, we combine the observations to quan- 
tify attributes of the host and flare relevant to the possi- 
bility they may be produced by a supernova or variable 
AGN. In Section 21 we attempt to account for the to- 
tality of these observations with supernovae and variable 
AGN hypotheses, but find that no known phenomenon 
other than tidal disruption is compatible with all of the 
observations. 

3.1. SDSS observations 

Figure [TU] shows the u, g, and r-band light curves for 
the flares (i.e., the difference images). Also plotted are 
the FUV and NUV fluxes from GALEX (see Sections 
13.21 fc 15751 for details on these GALEX observations), and 
the flux estimated from the CRTS optical observation of 
TDE2 3 months prior to the first SDSS observation in the 
flaring state (see Section l3.3.3l) . Both flares were detected 
by SDSS in the first observation of a Stripe 82 season, so 
were most likely past their peaks when first detected; this 
is confirmed by the CRTS detection in the case of TDE2 
(see Section f3. 3. 3[) . For comparison, we show fits of the 
SDSS data to f v (t) oc (t — toY ', where to is the time of 
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Fig. 10.— UV and optical light curves for TDEl and TDE2 as a 
function of days since observed peak. The SDSS difference image 
flux (i.e., host subtracted) in the r,g, and u bands is shown with red, 
green, and blue solid circles. The orange open square indicates the 
mean of the 3 CSS obser vation s 3 months before the first SDSS flare 
observation (see Section l3.3,31 , The dashed (solid) lines display the 
result of fitting a (t — t]j) p power law decay with p = —5/3 (—5/9) 
to the SDSS observations only. The corresponding NUV curves 
(purple) are obtained from the black body fit to the mean optical 
colors of the flare, assuming no cooling. Because the UV baseline 
of the host of TDEl is unknown, we show the GALEX aperture 
flux for TDEl, while for TDE2 we show the GALEX difference 
flux. 
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Fig. 11— The SED for TDE1 (top) and TDE2 (bottom). The 
optical baseline flux of the host is sho wn (black squares) with 
the best-fit combination of eigen spectra lIBlanton Ik. Roweial2007l i 
(grey line). For TDE2, the GALEX observations before the flare 
(black circles) are also used in this fit. The flux of the flare in 
the optical difference images is shown with red stars. It is well-fit 
by a black body shown as a blue line, with dotted lines indicating 
the 1-a uncertainty. The post-flare UV flux is shown with purple 
triangles; for TDE1 it is the total flux (no baseline being available) 
while for TDE2 the difference flux is shown. The detailed interpre- 
tation of the post-flare UV detections are unclear; the dynamics 
of the tidal debris and associated accretion are complex - a single 
black body need not be valid and the cooling behavior is uncertain. 



disruption. The SDSS optical observations do not cover 
a long enough period to break the degeneracy betwe en t n 
and p . For p = —5/3 (the fallback rate of the debris (jReesi 
1 1988T ) ) , the inferred time delay between when the disrup- 
tion occurred and when it was first observed by SDSS 
is 107 and 220 days for T DE1,2 respectively, while for 
p = -5/9, as predicted by iStrubbe fc Quataertl (|2009ft 
for the initial super-Eddington outflow phase, these be- 



come 24 and 55 days. In Fig. [TO] we show the p = —5/3 
and —5/9 extrapolations assuming constant temperature 
to obtain the NUV magnitude, for orientation. 

The photometric observations are well-fit by a rest- 
frame black-body spectrum with temperatures given in 
Table [U as can be seen in Fig. [TTJ An estimate of the 
cooling rate was obtained by least-squares fitting for the 
slope of color as a function of time, using only SDSS 
observations starting with the peak of the flare. Wc 
also computed the mean colors by averaging all obser- 
vations of the flare. Comparing the mean color to the 
cooling (Fig. [TOk). gives strong evidence that the TDE 
flares are not (ordinary) supcrnovae: they are very much 
bluer than any observed off-center flare in our sample and 
show negligible cooling, whereas SNc arc cither blue and 
rapidly cooling or red with little cooling because they 
start hot and cool very rapidly. Furthermore, Fig. ITOb 
shows that the colors of the TDE 1,2 flares are strikingly 
different from those of QSOs; if these were AGN flares, 
the spectrum of the flare itself (with the galaxy sub- 
tracted) should more closely resemble the spectrum of 
a QSO flare. 

To estimate the black hole mass (Mb h) of the host of 
TDE 1,2 we use the lHaring fc Rixj (|2004T ) black hole mass- 
bulge mass relation , calibrated for the SDSS r-band by 
iTundo et all (|2007l) . Unfortunately, the S/N in the co- 
added images of the host of the flare is too low to mea- 
sure the bulge magnitude by decomposing the images 
into a bulge and disk component. We therefore esti- 
mate the bulge magnitude using two different assump- 
tions for the bulge-to-total ratio (B/ T). (i) The typi- 
cal ra tio for SO galaxies, B/T = 0.55 (jAller fc Richstong 
12002} ) . (ii) The ratio that follows from the correlation be- 
tween B/T and the concentration index (i.e., R90/R50, 
where R90 and R50 are th e radii enclosi ng 90 and 50 
per cent of the galaxy flux) (IGadotti||2009f ). which yields 
B/T = 0.16,0.13 for TDE1,2. Using these estimates 
to obtain a range for the bulge luminosity of the host 
we obtain Mbh = (6 - 20) x 10 6±0 3 M Q for TDE1 and 
Mbh = (2 - 10) x 10 7±0 - 3 M Q for TDE2, where the error 
in the exponent reflects the scatter in the MeH-bulge lu- 
minosity relation and the range in the prefactor reflects 
the uncertainty in the bulge luminosity. 

3.2. TDE1 
3.2.1. Additional observations of TDE1 

For TDE1 we have a single GALEX NUV and FUV 
detection on 2008 October 25, ~ 800 days after the 
optical flare. The GALEX UV flux values, corrected 
for the energy lost due to the 6" radius aperture, and 
using the pipeline-generated sky background value in 
counts per pixel and assuming Poisson errors aperture 
(IMorrissev et al.|[2007h . are: FUV = 23.3 ± 0.3, NUV = 
23.0 ± 0.2. These UV fluxes are plotted in Figs. [TO] and 
1111 without host subtraction since we have no pre-flare 
observation to establish the baseline. 

We observed the host galaxy of TDE1 on 2009 Novem- 
ber 7 with the MagE spectrograph on the Magellan II 
Clay telescope with a spectral resolution of R = 4100 
or 1.6 A. The spectrum, shown in Fig. [TO] shows 
Balmer absorption features, as well as Mg lb and Nal 
absorption, which yields z = 0.136 ± 0.001. No emis- 
sion lines arc detected: Lrq < 4.4 x 10 37 ergs _1 and 
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Fig. 12. — a) The cooling time measured by fitting u-rasa function of time. The two TDE candidates are incompatible with being 
ordinary SNe because they show no signs of cooling and have very high u to g flux ratios, b) Color-color diagram usi ng the mean colors o f 
the decaying part of the light curve. Contours containing 90 and 95% of 14,776 nearby (z < 0.8) spectroscopic QSOs (Richards ct al. 2004) 
are also shown. Fewer objects appear in the upper panel because two simultaneous detections in both bands are required to measure the 
cooling time. For SNe and potential TDEs, the flux shown is that of the difference image; AGNs are shown in their high state. Colors are 
obtained from the error-weighted mean of all observations of the flare. Blue boxes mark flares from hosts that are not identified as AGNs 
based on their spectra or color, but whose variability in other seasons shows they are, in fact, AGNs. The supernovae in this work are 
selected pur ely geo metrically by being off-center (Eq. [SJ) and thus their properties are unbiased. The SNe that survive the TDE quality 
cuts (section l2. 2.31 1 arc indicated with a downward pointing triangles. 



TABLE 3 

Properties of the hosts of the TDEs. 



Flare 


SDSS ID 


RA 


Dccl. 


M r 


u-g 


g-r 










(J2000) 


(J2000) 










(Mq) 


TDE1 


J234201.40+010629.2 


350.95257 


-1.1361928 


-19.85 ±0.02 


1.95 ±0.3 


0.73 ±0.02 


(6 


- 20) x 10 6±0 - 3 


TDE2 


J232348.61-010810.3 


355.50586 


1.1081316 


-21.30 ±0.02 


0.99 ±0.2 


0.73 ± 0.05 


(2 


- 10) x 10 7±0 - 3 



Note. — The host magnitude and colors are obtained from the K-corrcctcd (Blanton & Rowcis 2007), inverse- variance- 
weighted mean Petrosian magnitude of the non-flare se asons. The black hole mass is estimated using the correlation between 
M r and A/bh (Harms & Rix 2004; Tundo ct al. 2007), using two different estimates for the bulge magnitude (see Section 

ED. 
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Fig. 13. — Host s pectrum of TDE1. A template 
Bruzual & Chariot (2003) spectrum is shown in red. Note lack of 
[OIII] emission lines and the Ha absorption. Hatch marks indicate 
the locations of the strong O2 telluric B- and A-band absorption 
features that were not removed in the data reduction. 

i[oni] < 1.3 x lO^ergs- 1 . 

3.2.2. Origin of UV emission of TDE1 

Star formation contributes little to the observed UV 
in TDE1. This follows from the Lho^fuv relation 
(jKennicuttl 119981 ) which predicts that if the FUV lumi- 
nosity observed by GALEX ^800 days after the opti- 
cal flare were due to star formation, the Ha luminos- 
ity should be 2.4 x 10 40 erg s _1 - two orders of mag- 
nitude above the observed upper limit. The contribu- 
tion of stellar sources that evade the correlation between 
FUV and Ha (e.g., blue horizontal branch stars) is con- 
str ained by using the color-magnitude relation observed 
by lHaines et al.l (|2008[ ) for galaxies that have been spec- 
troscopically identified as passively evolving. At the lu- 
minosity of the host of TDE1 (M r = -19.85), these 
galaxies are observed at NUV — r = 5.4. The scatter 
in this relation, cnuv-i- = 0.37, can be attributed to dif- 
ferent amounts of residual star formation. Thus the UV 
detection of the host of TDE1, with NUV - r = 3.6, 
is a 4.6-cr blue outlier to these galaxies. Moreover, 
with FUV — r = 3.7±0.3, it is distinctly bluer than 
passive galaxies of similar luminosity, which cluster at 
FUV - r = 7. The UV color, FUV - NUV = 0.0 ± 0.36, 
is also bluer than any of the early type ga laxies which 
have been targeted by the SAURON project dJeong et al.l 
l2009t) to study residual star formation. These galaxies 
are observed to cluster at FUV — NUV = 1.5, with the 
bluest object at FUV — NUV = 0.6. Recent or residual 
star formation is thus excluded as a significant contribu- 
tor to the observed UV emission detected 800 days after 
the first observation of TDE1. 

Nor can an active nucleus account for the observed 
UV flux 800 days after the TDE1 flare. The upper 
limit on the [OIII] luminosity from the Magellan spec- 
trum places a limit on the optical baseline luminosity of 
a possible active nucleus. Using the conversion between 
[O III] and optical lumin osity for type 1 AGNs obtained 
bv lHeckman et aL1(|2004| ). yields an upper limit on the lu- 
minosity of an active nucleus in the TDE1 host: £5000 < 
4.2 x 10 40±a34 ergs -1 , where -L5000 is the monochromatic 
continuum luminosity at 5000 A i n the rest frame. Al- 
though the iHeckman et al.l (|2004l ) relation is not valid 



for [OIII] emission from LINERs (IHeckman! I"l980l ). we 
can still use it to obtain an upper limit on the baseline 
luminosity of a potential AGN in such galaxies, because 
at similar [OIII] luminosity, the active nuclei of LINER 
galaxies have an order of magnitude lower bolometric lu- 
minosity compared to normal AGNs, as well as an lo wer 
optical lum inosity for a given bolometric luminosity (|Hol 
1999, 2004). We convert £5000 to a g-band magnitude 
(centered at A = 4670 A), by assuming that the luminos- 
ity per unit wavelength at the g-band is equal to £5000- 
This is the conservati ye approach since F„ oc t /~ 44 is 
typical for an AGN (jVanden Berk et al.l 12001 . This 
yields NUV - g < -2.9 ± 0.9, a NUV to p-band ratio 
greater than the maximal value that can be reached with 
photons that are in thermal equilibrium; such an extreme 
color has never been observed for a ny of the sources de - 
tected in both SDSS and GALEX (jBianchi et al.l 120071 ). 
Because we obtained our upper limit on the g-band lumi- 
nosity of the accretion disk from the [OIII] line, this limit 
applies to the baseline luminosity of the AGN, i.e., time- 
averaged over the light crossing time of the nar row line 
region, > 10 2 yr (jMuravama fc Taniguchil fl998h . Thus 
we conclude that the UV flux present 800 days after the 
optical flare does not originate from the baseline of an 
AGN. 

3.3. Additional observations of TDE2 

We have GALEX NUV and FUV detections of the 
TDE2 host in both the pre-flare (2003 24 August) and 
post-flare (2008 14 October) state. The TDE2 host was 
also observed by GALEX on 2008 30 August, but the 
image was at the edge of the field and the sky was noisy, 
so we display the NUV value in Fig. ITOb ) but do not use 
it further. Figure [14] shows the GALEX images. The 
pre-flare GALEX UV flux value s, corrected for the en - 
ergy lost due to the 6" radius ([Morrissev et al.l [20071 ). 
are FUV = 23.0 ± 0.3, NUV = 22.8 ± 0.2. The GALEX 
post-flare difference magnitudes are: FUV = 21.4 ±0.2, 
NUV = 21.1 ± 0.1. Fig. [TTJd shows the optical-UV SED. 

In addition to the UV data, we have two serendipitous 
and two follow-up optical observations of TDE2 and its 
host as discussed in greater detail below. The spectrum 
of the host galaxy was obtained in a follow-up obser- 
vation with the William Hcrschel Telescope (WHT) in 
November 2010. Serendipitously, the first points on the 
SDSS light-curve of TDE2 suggested it might prove to b e 
a supernova and an alert was issued (|Poimans 12001 . 
As a result, an optical spectrum was obtained with ESO's 
New Technology Telescope (NTT) on 2007 September 
18, 4 days after the first detection of the flare, and radio 
observations were obtained with the VLA 7 and 92 days 
after the first SDSS detection; these observations are dis- 
cussed in Section 13.3.11 below. Furthermore, a review of 
Catalina Real-time Transient Survey observations after 
discovery of TDE2 reveals a detection of the flare 95 days 
before the beginning of the Stripe 82 observing season. A 
preliminary report on this data is given in Section 13.3.31 
a more thorough analysis using also a Keck spectrum of 
the host galaxy obtained in October 2010 will be reported 
elsewhere. 

3.3.1. TDE2 Host and Flare Spectra 

In November 2010, 850 days after the first SDSS detec- 
tion of the flare, the host galaxy of TDE2 was observed 



Optical discovery of tidal disruption flares 



13 




Fig. 14.— GALEX FUV and NUV images before and during 
the flare in TDE2. Red circles show the 6 arcsec radius aperture 
used to measure the photometry. The plate-scale of the image is 
1.5 arcsec/pixel. Note that in the second epoch (2008-08-30), the 
source is at the edge of the field where the photometry is more 
susceptible to systematic errors due to distortions in the PSF. 



with the WHT for two nights. On 01 November two 
ACAM exposures of 900s were taken using a V400 grat- 
ing and a slit width of 1". On 05 November, the host 
was observed again using the same specifications and ex- 
posure times, except that a slit width of 1'.'5 was used. 
The data of both nights (lh total integration time) was 
reduced separately, using standard IRAJ0 routines, to 
yield four ID spectra which were then combined. The 
overall flux normalization was obtained using a standard 
reference star; the wavelength dependent flux normaliza- 
tion was fixed using the SDSS PSF flux of the host. 

The NTT spectrum - taken 5 days after the first 
SDSS detection of the flare - wa s taken using the E SQ 
Multi-Mode Instrument (EMMI) (|Dekker et al.lll986D in 
the RILD mode (Red Imaging and Low-Dispersion spec- 
troscopy). The grism that was used, Gr 2, has a wave- 
length coverage of 3800-9200 A, 300 grooves/mm, a wave- 
length dispersion of 1.74 A/pixel, and a spatial resolu- 
tion of 0'.'166 /pixel (before binning). During the obser- 
vation a binning of 2x2 was used. The slit width was 
1'.'5 and the exposure time was 1200 s. The flux cal- 
ibration was performed using spectra of spectrophoto- 
metric standard stars observed with a 5" slit. The stan- 
dard star spectra were also used to construct a telluric 
spectrum by isolating the telluric absorption in them. 
The spectrum of TDE2 was then corrected by subtract- 
ing a scaled version of the telluric spectrum. Finally 
the spectru m was corrected fo r dust ext inction in the 
Milky Way (ISchleeel et al.lH998h using the lCardelli et al.l 
(|1989t) CCM extinction law with R v = 3.1. For a more 
detailed description of the data reduction we refer to 
IQstman et al.l (|2011l) . Because the spectrum of TDE2 
was not obtained in parallactic angle, it is affected by 
differential atmospheric refraction. The wavelength de- 
pendent flux loss was estimated and corrected lor. tak- 
ing into account the seeing, slit width, airmass, the angle 
from parallactic angle and the wavelength a t which the 
object was centered on the slit (|Owenslll967[) . 

Figure [TS] shows the host and flaring-state spectra to- 
gether with the SDSS photometry; a zoom on Ha is 

13 IRAF is distributed by NOAO, which is operated by AURA, 
under cooperative agreement with NSF 



shown in Fig. 1161 The H a line luminosity in the host 
spectrum, measured by fitting the observed flux to a 
Gaussian whose width is given by the PSF, is Lhq = 
1.2 ± 0.2 x 10 40 ergs _1 ; we simultaneously fit for the 
[Nil] emission lines which yields log([NII] A6583/Ha) = 
— 0.4±0.2. Using the redshift obtained from this simulta- 
neous fit as an initial guess, we measure the flux of other 
well-known emission and absorption lines in the spec- 
trum; only lines measured above the 3-tr level, that are 
displaced by no more than 2% from the initial redshift, 
are considered real detections. The results of running 
this procedure for the host and flare spectrum are shown 
in TableH Using all detected lines in the NTT and WHT 
spectrum we obtain z = 0.2515±0.0036. From the WHT 
spectrum we obtain i[oinj < 3.6 x 10 39 ergs _1 . 

Because TDE2 occurred at the center of its host, the 
NTT spectrum contains a large galaxy component that 
has to be subtracted to obtain the "pure" flare spec- 
trum. After subtraction, no narrow (i.e., unresolved) 
lines remain (Fig. [16l bottom panel). However at the 
rest-frame wavelength of H Q we identify a broad feature 
which can be fitted by a single Gaussian of a = 75 ± 5 A 
or FWHM = 8 x 10 3 kms~ 1 and has an equivalent width 
of 87 ± 5 A. 

Besides the NTT spectrum, the initial identification of 
TDE2 as a possible SN also triggered radio observations. 
We reduced the data from two 8.5 GHz VLA observa- 
tions, obtained as part of project AS 887 - the only VLA 
observations of this field after this event. The first VLA 
observation was 7 days after the first SDSS detection of 
TDE2 and the second was 85 days later. For both ob- 
servations, we calibrated the recorded visibilities to flux 
density using data from short observations of 3C48, and 
used visibility data of P MN J2323-0317 for phase calibra- 
tion (jBaars et al.l[l977 ). Imaging the two observations 
separately yielded no detection of TDE2 at either epoch. 
Combining the visibility data of TDE2, we obtained an 
image with an RMS noise of ~ 35 /xJy per beam using 
natural weighting. No source was detected at the loca- 
tion of TDE2, allowing us to place a 3-cr upper limit of 
^0.1 mJy on its 8.5 GHz flux density. 

3.3.2. Interpretation of Host and Flare Spectra 

The WHT host spectrum is consistent with the conclu- 
sion based on photometry in Section l2.3.3[ that the flare 
is not due to extreme variability of a hidden AGN. Such 
an AGN would reveal itself by a higher [Nil] A6583 to Ha 
ratio than observed. Using the BPT diagram (Eq. [3]), we 
see that the observed ratio is consistent with the narrow 
Ha origi nating entirely fr om star formation. Moreover, 
using the lKennicuttl (|1998[ ) relation, we conclude that the 
star formation rate inferred from the narrow Ha lumi- 
nosity (0.5 Mq yr _1 ) is consistent with the rate expected 
from the pre-flare FUV luminosity (0.7M Q yr^ 1 ). 

Employing the wid ely-used cross correlation method 
([Tonrv fc Davis! fP379l) fo r identifying SNe, using SNID 
( Blondin fc Tonrvl 120071 ) on the host-subtracted spec- 
trum, gives the best match to be an early SN type Iln 
spectrum, but this match is not very convincing: the de- 
gree of cross-correlation of the best-matching SN spec- 
trum is similar to the value obtained for a template 
galaxy or AGN spectrum and is close to the cut-off value 
of the SNID software. A further comparison of TDE2 to 
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TABLE 4 
Narrow Lines 



Line name EW host EW flare Flux host Flux flare 

(A) (A) (lCr^ergs^cm- 2 ) (lCT 16 crgs^cm- 2 ) 



K a 4.2 ±0.8 4.1 ±0.7 0.8 ±0.1 1.8 ± 0.3 

[Nil] 6585 1.4 ±0.7 2.8 ± 0.7 0.2 ±0.1 1.3 ± 0.3 

Ca H -5.7 ±1.8 -2.8 ± 0.7 -0.8 ± 0.3 -1.5 ±0.4 

Ca K -7.6 ±1.3 -3.3 ±0.7 -1.1 ±0.2 -1.9 ±0.4 

H { -3.4 ±1.1 -0.6 ±0.2 



Note. — Narrow (i.e., unresolved) lines detected in the host (WHT) and flaring- 
state (NTT) spectrum. Line fluxes are obtained by fitting a Gaussian profiles whose 
width is given by the PSF to the observed flux. Only lines detected above 3cr are 
considered. 



Flaring-state spectrum (NTT), 2007-09-18 




4000 5000 6000 7000 

Rest — frame wavelength (A) 

Fig. 15. — Host and flaring-state spectra for TDE2: NTT spec- 
trum in the flaring state (black) and WHT spectrum of the host 
(red). We also show the SDSS PSF UStoughton et al.|[2002T) flux 
of the host in the quiescent state (open squares) and the SDSS 
flux during the flare (purple stars). We see that the photomet- 
ric and spectroscopic flux calibration agree reasonably well. The 
narrow Ha emission is consistent with the level expected from star- 
forma tion as implied by the pre-flare FUV luminosity l IKennicutH 
1199811 . A zoom on Ha is shown in Fig. 1161 

type Iln SN is presented in Section |4"T2"1 

An exciting signature of a TDE is line emission from 
tidal disruption debris which is shifted by the very high 
Kcplcrian velocities in the vicinity of the supermassive 
black hole. However, the strength of this line emission 
is predicted to be much weaker than the flare continuum 
iBogdanovic et ail ()2004l ). and is likely below the detec- 
tion threshold of the TDE2 flare spectrum. We note that 
the narrow line at A A ss 30 A (1.5 x 10 3 kms" 1 ) blueward 
of Ha (Fig. [TBI top panel) is a candidate to be such a 
shifted line, but it is detected only at the 2cr level, with 
an equivalent width of 1.5 ± 0.7 A. 

The interpretation of the intermediate width line in the 
host-subtracted flare sp ectrum (Fig. [TBI bot t om p anel) 
is not straightforward. iStrubbe fc Quataertl (|2011[ ) cal- 
culate that the line-broadening due to the line-of-sight 
velocities of the super-Eddington outflows is probably 
too broad to be detectable. However the stellar debris 
orbiting the black hole at e ccentric orbits could a lso pro- 
duce a broad component l)Komossa et al.l I2009D . Since 
intermediate width H Q lines arc a well-known property 
of many AGN and some type II SNe, their detection in 
the host-subtracted spectrum of TDE2 is not particularly 




6400 6450 6500 6550 6600 6650 6700 



Rest-frame wavelength (A) 



Velocity (kms l ) 

, , -20000 -10000 10000 20000 

1.6 p — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 




6200 6400 6600 6800 7000 

Rest-frame wavelength (A) 



Fig. 16. — Zoom on the Ha emission of TDE2 for the flaring- 
state and host spectrum (top) and the host-subtracted spectrum 
(bottom). The host spectrum is arbitrarily rescaled for compari- 
son to the shape of the flare spectrum. For display purposes, all 
spectra are smoothed using a Gaussian filter with a = 1.5 A. The 
intermediate width Ha line in the host-subtracted spectra can be 
fit with a single Gaussian to find a FWHM of 8 X lO^kms^ 1 . The 
feature at 6529.7 A is detected at the 2a level (EW = 1.5 ± 0.7 A); 
it may be interpreted as weak, blueshifted Ha emission (with a 
velocity of 1.5 X 10 3 kms -1 ), but could also be an artifact. The 
full range of the host-subtracted spectrum is shown in Fig. 1211 

constraining to the nature of this flare. 

3.3.3. CRTS observations 

The Catalina Real-Time Transient Survey (CRTS) an- 
alyzes data from the Catalina Sky Survey which repeat- 
edly covers >30,000 square degrees on the sky, in order 
to search for optical transients with timcscalcs of minutes 
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to years dDrake et all 120091) . The Catalina Sky burvey 
(CSS) Schmidt Telescope and Mount Lemmon (MLS) 
telescopes are located north of Tucson, Arizona and sur- 
vey the northern sky using unfiltcred 4k x 4k CCD cam- 
eras. Images from the MLS telescope cover 1.1 deg 2 and 
reach V =21.5, while images from the CSS cover 8 deg 2 
and reach V = 19.5. On a clear night, these two tele- 
scopes cover ~ 1500 deg 2 of sky in sequences of four 30 
s exposures. Although CRTS began in 2007 November 
archival CSS and MLS data dates back to 2004. 

CRTS observations of TDE2 are available both in be- 
tween and beyond the SDSS observational seasons. This 
data and a Keck spectrum of the host will be presented 
in a forthcoming paper, with full calibration and opti- 
mized difference imaging so they can be quantitatively 
combined with other observations. Here we only show 
the mean "unfiltcred" CSS difference magnitude (within 
~ 20% of the SDSS r-band) of the three CSS detections 
of the TDE2, obtained 95 to 80 days prior to the first 
SDSS detection. These data are shown in Figure [TOl to: 

• date the first SDSS observation and NTT spectrum 
as being at least 95 and 100 days respectively after 
the peak of the flare, and 

• show that (Fig. ITOj) the peak luminosity of the flare 
is considerably higher than observed by SDSS, by 
approximately one magnitude. 



4. COMPARISON OF TIDAL FLARES TO AGN 
FLARES AND SUPERNOVAE 

In this section we quantify the likelihood that TDE1 
and TDE2 are examples of some already-observed phe- 
nomenon. 

4.1. Comparison of TDEs to AGN flares 

Our pipeline was designed to e xclude AG N f lare s based 
on the host properties (Sections 12.31 13.21 and I3.3|) . Here 
we examine the properties of the flare observations them- 
selves to understand how often the behavior displayed by 
TDE1,2 may occur by chance in the AGN population. 
We use two attributes of variable AGNs that we mea- 
sure in our data - the range of variability per season and 
the spectrum of the most extreme flares over the entire 
Stripe 82 observing period - to quantitatively compare 
TDE1,2 to AGNs. First, we consider the likelihood for 
an AGN to have a flare meeting our selection require- 
ments in one season, yet be as quiet as TDE1,2 in the 
other observed seasons. Next, we consider the spectrum 
of flare amplitudes in the AGN population and compare 
to those of TDE 1,2. Finally we compare the spectral 
properties of the TDE1,2 flares to AGN flares. 

4.1.1. AGN Variability 

Probability of an isolated flare episode in QSOs — As dis- 
cussed Section |2. 3. 3[ we rejected three candidates which 
pass all cuts except that their flux in other seasons is not 
consistent with being constant . We can use the measure 
of variability developed there to determine the likelihood 
that an AGN is as quiet as TDE1,2 in all but one ob- 
serving season. We have a sample of variable AGNs with 
the same selection criteria as TDE1,2 - our sample of 
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Fig. 17. — The x 2 /DOF distribution of all non-flare seasons 
with three or more observations (x?)> for the spectroscopic QSOs 
and Seyfert galaxies. The values of Xs m the non-flare seasons of 
TDE1.2 are respectively {1.4, 2.6, 0.6, 1.7, 1.9} and {1.0, 1.6, 1.0, 
2.0}. The inset zooms in on the first few bins of the histogram, 
using a linear bin size. 



flares that are in QSOs and identified Seyfert galaxies - 
with which we can quantify the flux excursions in "off" 
seasons. To do this, we introduce x 2 '- the x 2 /DOF per 
season, s. The median flux excluding the season that 
contains the main flare is used as a model for the light 
curve in other seasons. We calculate xl f° r an seasons 
(other than the one with the primary flare) having three 
or more detections, since that is the minimum number 
of detections for the TDE candidate sample. Fig. IT71 
shows the distribution of these x 2 values. Let f^' be 

the fraction of AGNs with as low or lower value of x 2 as 
TDE1,2 in season j. Our estimate for the likelihood for 
an AGN to have as little activity in the off seasons as 
displayed by TDE 1,2 is then 



N, 



^<II 2 /: 



(1,2) 



(4) 



Here the product runs over the non-flare seasons with 
three or more detections; N s = 5,4 for TDE 1,2. The 
factor 2 is inserted in Eq. (j4} because the mean value 
of fj for AGNs is 0.5. Using the relative photometry 
light curves introduced in Section 12 .31 we find Pagn _ 
2 x 10~ 6 , 2 x 10~ 5 for TDE1.2 respectively. 

This estimate of the chance probability for a variable 
AGN to have several quiet years surrounding a major 
flare, and thus to be able to mimic TDEl,2's variabil- 
ity properties, assumes the flux variability in the years 
surrounding the major flare is uncorrelated. This is a 
reasonable first approximation and can be improved by 
studying the statistics of the AGN fluctuations in the 
years near a major flare. It is important to emphasize 
that we are not making the assumption that the natural 
flux variability of AGNs in seasons near a major flare is 
the same as in randomly chosen seasons long before or 
after a major flare: if there is an enhancement or suppres- 
sion of flux excursions before or after major AGN flares, 
it is captured in our measured x 2 distribution, which is 
measured in AGNs in years surrounding a flare selected 
by the same criteria as for TDE 1,2. 
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Fig. 18. — The flux increase of QSO flares with respect to the 
baseline of the flare for the g, r and i-bands. We fit the histogram 
with P(AF/F) oc (AF/F) a for all bins with a flux increase larger 
than 10% and obtain a = —3.9 ± 0.2. The probability of finding a 
flare as large or larger than a given AF/F is obtained by integrating 
this fit and multiplying by 112/1304, the ratio of flaring QSO to 
all QSO that have been searched for flares. For reference, we also 
give the lower limits on the flux increase of TDE1,2 with respect to 
a hypothetical AGN baseline flux as derived from the (upper limit 
on) the [OIII] luminosity. 

Probability of comparably large flares in QSOs — We can 
use another property of the flaring QSO sample to obtain 
a second, independent probability measure that TDE1 or 
TDE2 are AGN flares. For this, we quantify the spec- 
trum of flux increase at the peak of QSO flares rela- 
tive to their baseline flux. An upper limit on the [OIII] 
line luminosity of a galaxy can be converted to an up- 
per limit on the baseline luminosity of its active nucleus 
at 5000 A in the rest-frame (jHeckman et al.ll2004ft . The 
observed luminosities of TDE1,2 at this wavelength im- 
ply a minimum flux increase with respect to the baseline 
state of the accretion disk of a factor 87, 21 respectively. 
A flux increase of this magnitude is extremely unlikely 
for an AGN: out of the 1304 extended QSOs in Stripe 
82 that we monitored for flares, the largest flux increase 
measured in the g, r or i bands, is a factor of 5. The 
spectrum of AF/F for the QSO flares meeting our se- 
lection criteria is shown in Fig. 1181 it is a power-law 
with slope a = —3.9 ±0.2. Using this power-law fit 
to calculate the probability of as large a flux excursion 
as seen in TDE1,2, if they were variable QSOs, gives 
P(AF/F > 87) = 3 x 10~ 7 , P(AF/F > 21) = 2 x 10~ 5 , 
respectively. 

Applicability of the above estimates — The probability es- 
timates above are valid if there is only one type of AGN 
variability. Evidence that a single phenomenon is re- 
sponsible for the variability observed in most accreting 
super-massive black holes can be found in the litera- 
ture. Usin g ~ 9000 spectroscopi cally confirmed QSOs of 
Stripe 82, iMacLeod et al.l (|2010ft showed that a damped 
random walk model can explain quasar light curves at 



an impressive fidelity level (0.01 — 0.02 mag), indicating 
that one single physical process, e.g., turbulent m agnetic 
fields within the accretion disk dKellv et al. 2009), is the 
dominate source of the variability. This analysis does 
not apply to BL Lacs, which show larger amplitude fluc- 
tuations (up to A m ~ 1) than QSOs on all timescales 
l|Bauer et al.|[200l , but as we shall see, the properties of 
TDE1,2 do not place them in this class. 

The violent fluctuations in BL Lac type AGNs are 
very likely due to a fluctuating contribution from a rel- 
ativistic jet. A g o od ex ample is the flare found by 
iVanden Berk et all (pOOl . The first SDSS photomet- 
rical observations showed a red source (g — r = 0.3); the 
SDSS spectrum obtained about a year later, showed this 
sources had faded by ~ 2.5 mag and re vealed a galaxy 
spectr um. Followup observations by IGal-Yam et al.l 
(|2002|) showed broad H Q in the host spectrum, and a 
second flare with a blue continuum. The radically dif- 
ferent SED for the second flar e, combined with a de- 
tection at 1.4 GHz in FIRST (jBecker et al.l 119951) and 
othe r radio catalogs that r evealed radio flux variability, 
lead IGal-Yam et al.l (|2002f ) to conclude the flare origi- 
nated from a radio-loud AGN, probably in the B L Lac 
class. We note that the IVanden Berk et al~l (|2002f ) flare 
is somewhat similar to Rl and R2, the flares we rejected 
based on additional variability of the host fsection r2.3.3l) . 
which are also detected in FIRST at i<i.4GHz > 90 mJy 
and have similar re d color s. TDE1,2 are nothing like the 
IVanden Berk et aLl ([2002D flare: their post-flare spectra 
show no broad lines, the flare SED is nearly constant, ad- 
ditional seasons of observations in Stripe 82 show no ad- 
ditional variability, and they are not detected in FIRST 
(or targeted VLA observations for TDE2). 

To quantify the differe nce between the host of the 
IVanden Berk et al.l (pOOl flare and TDE1,2, we use 
the 3-ct upper limit on the radio luminosities (L„ < 
1 x lO^egss^Hz- 1 at 1.4 GHz for TDE1 from FIRST 
and L u < 2 x 10 29 egs s^Hz -1 at 8.5 GHz for TDE2) 
to compute the optical-to-radio spectral indices, a ro = 
-log(F r /F )/log(j/ r /V ) < 0.08,0.04 for TDE1,2 using 
the g-band peak flux for F Q . From SDSS observations, 
optically selected radio-loud BL Lacs ha ve a Gaussian 
distrib ution with (a ro ) = 0.42 ± 0.08 (jPlotkin et all 
l2010bf) . hence the low radio-to-optical ratio of the hosts 
of TDE1 or TDE2 is not consistent with a blazar-origin 
of these flares; they would have to members of a new 
class of radio-quiet BL Lacs with violent flares. 

About 80 radio-quiet objects (a TO < 0.2) with spectra 
that resemble BL Lacs at optical wavelengths (i.e., no or 
very wea k emission lines, blue continuum) , are known at 
z < 2.2 (IPlotkin et al. l l2010allbl ). Using Stripe 82 data, 
IPlotkin et all (|2010aD conclude that the level of optical 
variability of these AGNs is consistent with other radio- 
quiet quasars: x 2 /DOF > 10 for all Stripe 82 obser- 
vations in the g-band (which can be compared to Fig. 
[8]) and no excursions from the mean flux greater than 
Am = 0.5. Hence all of these peculiar AGNs would be 
identified by our pipeline and excluded based on their ad- 
ditional variability, and none show flares that are large 
enough to hide the underlying AGN. We conclude that 
on the radio-quiet branch of AGN, there is no evidence 
for a different, more violent mode of variability as seen 
in radio-loud BL Lacs, but more data is needed to fully 
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Fig. 19. — Colors of spectroscopic QSO. The quiescent state is 
calculated from the mean PSF flux of the non-flare nights. The 
high state is obtained by adding the flux of the quiescent state to 
the difference image and taking the mean colors. No color change 
during a flare is indicated by the striped blue line. We see that 
QSOs in their high state are slightly bluer; this effect is most pro- 
nounced for flares from relatively red QSO (quiescent u — r > 0.5) 
which can be explained by stellar contamination. For reference, 
the u — r colors of TDE1,2 measured in the difference image are 
-0.62 ±0.06, -0.36 ±0.03. 



exclude the existence of such a mode. 

4.1.2. Colors of TDE and QSO flares compared 

Figure IT2b shows that the colors of the TDE flares, 
as measured in the difference image, fall outside the lo- 
cus that contains 97% of all spectroscop ically confirmed 
low-redshift (z < 0.8) QSOs in SDSS (ISchneider et al.1 
[2001 . Here we investigate whether, when they flare, 
QSOs may change color radically enough that our TDEs' 
SEDs could be consistent with that of a flaring QSO. 
To explore this question, we compare in Fig. Qj)] the 
u — r colors of QSOs in the flaring and baseline states. 
QSOs are bluer in their high state than in the quiescent 
state, but the change is not nearly enough to push a QSO 
out of the QSO-locus sufficiently to match the colors of 
TDE1 or TDE2. A large color change is observed only 
for flares from QSOs with a relatively red baseline, i.e., 
where the galaxy contribution is relatively large and is 
not subtracted, whereas for the TDE flares the galactic 
contribution has been subtracted. 

4.1.3. Comparison of TDEs and Seyfert flares 

Comparing the properties of the TDE1,2 flares to the 
flares of QSOs is powerful because the large number of 
QSOs allows the dispersion in the QSO flare properties 
to be measured. However the TDE hosts are clearly 
not QSOs so we also compare specifically to the three 
examples in our flare sample in which the host galaxy 
is a Seyfert as defined in Section 12.3.11 (galaxies hav- 
ing broad lines or satisfying the BPT criteria, not la- 
beled QSO by SDSS). In Seyfert galaxies, the AGN con- 



tributes only a moderate fraction of the total flux, un- 
like in QSOs for which the AGN dominates by defini- 
tion. Hence for Seyferts the difference image flux gives 
a better measure of the flux of the AGN, and we there- 
fore compare the properties of the difference images of 
the TDEs with those of the three flares with identi- 
fied Seyfert hosts which pass our selection criteria. The 
magnitude and color decay rates of the three Seyfert 
flares cluster around zero; the bluest Seyfert flare has 
u — g = g — r = —0.2 which is roughly as blue as TDE2. 
To obtain more thorough comparison of the broadband 
properties, we fit each of the five mean SEDs (i.e., the 
average color) of the difference image to black body spec- 
tra. The Seyfert flares all have poor fits to a black body 
(X 2 /D0F = 5.9, 7.6. 22.7) while TDE1,2 are both well-fit 
by black-body spectra (% 2 /DOF = 1.9, 0.4), at optical 
wavelengths. 

The light curves of the Seyfert flares are also different 
from those of TDE1,2. They are more symmetric (i.e., 
rising and falling at similar rates) or show more substruc- 
ture, and they also show variability in the non-flare sea- 
son. To illustrate this, the light curve of a Seyfert galaxy 
that hosted at nuclear flare is shown in the bottom right 
panel of Fig. [7] 

4.2. Comparison of TDE1,2 to Supernovae 

In this section we show that TDE1 and TDE2 are un- 
like any off-center SN recovered by our pipeline, based on 
the properties of the flares. We examine the possibility 
that TDE1,2 could be exotic SNe, but find no examples 
of SNe which resemble the TDEs in all significant re- 
spects. We calculate the likelihood that TDE1,2 are the 
most similar type of known SNe - SNIIn - and are by 
chance as close as observed to the centers of their galaxies 
to be < 0.08 %. Finally, we consider the possibility that 
TDE 1,2 are examples of a new type of extreme SN which 
occur only near the centers of galaxies, but this requires 
a thousand-fold or greater enhancement in the rate of 
such explosions in stars in the nuclear region compared 
to elsewhere. 

4.2.1. Constraints from late-time UV emission 

The residual UV radiation detected 781, 485 days after 
the first optical detection of TDE1,2 is incompatible with 
observed UV properties of SNe. To estimate their max- 
imum possible NUV flux at these times, if TDE1,2 were 
SNe, we assume that when first detected in the optical, 
they had the bluest NU V to q-band ratio eve r observed in 
any SN: NUV-.g = -1 (|Bianchi et al.ll2007h . We further 
assume that the UV flux falls off with the slowest linear 
decay rate of UV magnitude that has been measured for 
a large sample of la, Ib/c and type II SN ([Brown et al.l 
120091 ): dTONuv/rft = 0.05 mag day" 1 , which is similar to 
the UV decay rat e of SN 2008es (jGezari et al.l l2009at 
iMiller et all 120091) . This yields upper limits on the 
NUV fluxes for TDE 1,2 under the SN hypothesis which 
are far below those observed: apparent magnitude of 
tonuv, SN > 59,42 compared to tonuv.tde = 23.0,21.1. 
We see that a linear decay, as commonly observed for 
normal SNe, vastly under predicts the UV flux, while a 
power law dec ay, which is predicted for stellar tidal dis- 
ruption events (Rees 1988), can explain the large residual 
UV flux, years after the flare event. 
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Fig. 20. — The change in magnitude and color measured after the 
peak of the flare. In this plot we only show the flares that passed 
the quality cut for the final TDE study (see Section [2.2.311 . The 
decay rate (color evolution) is obtained from the slope of the best- 
fit line to the flare magnitude (color) after the peak as function of 
time. All identified SNe show significant cooling, as expected, and 
their light curves decay faster than TDE1 and TDE2, except for 
one: a type II-P SN whose mean color, u — g = 1.3, is far redder 
than the TDEs. One flare from the photometric AGN sample is 
more like a SN than an AGN flare; this is consistent with the 0.9 
SNe expected in the sample for the TDE analysis. 



The sample of iBrown et al.l (|2009f) does not include 
type Iln SNe. To obtain an estimate of the limit 
of the NUV magnitude for this type of stellar explo- 
sion we use t he slo west S-band decay measured for the 
iKiewe et al.1 (|2010f ) sample of type Iln SNe (dms/df = 
0.016 magday" ), to find mNuv.iin > 32,26. This up- 
per limit on the NUV flux for type Iln SN is orders of 
magnitude below the observed late-time UV flux for both 
TDE1 and TDE2. 

4.2.2. Comparison of TDE1,2 to flares in our SN sample 



In Section [2 .2. 2 1 wc defined a sample of 85 flares which 
are clearly off-center from their host (Eq. [2]). This sam- 
ple allows up to compare the properties of TDE1,2 to 
normal SNe (exotic SNe are discussed below). Because 
this sample is obtained by a cut on the host-flare distance 
only, their properties are representative of the properties 
of SNe that appear in the nuclear sample, which would 
not be the case if we compared to SNe identified in spec- 
troscopic follow-up campaigns. We note that requiring 
at least three it-band detections, as is imposed for the 
potential TDE sample, biases the selection toward blucr- 
than-average flares, likely increasing the fraction of type 
Iln SNe. Even so, the TDE1 and TDE2 flares are sig- 
nificantly different from all SNe flares we observed with 
respect to the mean colors and cooling rate measured us- 
ing the mean u — g, g — r colors and the slope of u — r, 
respectively (Figs. IT2"j) . Also the decay rate measured 
using the slope of g is lower for TDE1,2 compared what 
is observed for the SN in the sample that meet the TDE 
quality cuts (Fig. |211 . except for one: a II-P SN (2006fg, 
iD'Andrea et al.ll2010D . However, as expected, this SN de- 
cays much faster in the w-band and is also red, u— g = 1.3. 

4.2.3. Comparison of TDEs to type Iln and Exotic SNe 



We saw above that the colors and cooling rates of the 
flares in our SN sample do not resemble those of TDE1,2, 
but this does not apply to all type Iln SNe and some 
exotic SNe. Here we review such SNe and compare them 
to TDE1,2. 

SNe Iln from the Caltec h Core-Collapse Project — Re- 
cently, IKiewe et al.l (|2010l) present ed four type Iln from 
the Caltech Core-Collapse Project (|Gal-Yam et al.ll2007T) 
which aims at producing an unbiased sample (i.e., es- 
sentially every young core-collapse SN that is observ- 
able was foll owed- up). At the ir peak, all four SNe 
presented by IKiewe et al.l (|2010D have similar or even 
bluer B — V colors and a similar decay rate as TDE 1,2. 
Two of them, 2005bx and 2005cp, show this decay rate 
over the full extent of their light curves (observed for 
50 and 119 days after their peak, respectively), while 
for the other two the flux drops rapidly after about 2 
months. SN 2005bx shows rapid cooling (d(B — V)/dt = 
0.04 day - ), while SN 2005cp has nearly constant color, 
d(B - V)/dt = 0.002 day -1 , compared to d{B - V)/dt = 
0.003, 0.004 day -1 for TDE1,2. The presence of one 
amon g four type Iln SNe in the small but unbiased sam- 
ple of IKiewe et al.1 (|2010T ) . which has similar colors and 
cooling with respect to the TDEs, thus shows that flare 
colors and cooling alone (Fig. [12"]) is insufficient to reject 
all SNe. 

We now compare the lumi nosity and spectra l fea- 
tures observed for TDE2 to the IKiewe et all (|2010D type 
Iln SNe which most resemble them. With Mb = 
-18.9,-18.0, SN 2005bx, cp are significantly less lu- 
minous at the peak than TDE2, whose equivalent Mb 
value is —20.3 when first observed, an unknown time af- 
ter the peak. In Fig. [21] we compare the host-subtracted 
flare spectrum of TDE2 to the spectra of SN2005bx and 
2005cp. The emission lines and P-Cygni absorption pro- 
files, a classic feature of SNe spectra, are absent in the 
TDE2 spectrum. While there is substantial dispersion 
within the sample of spectra of SN Iln, at all times in 
their evolution, the SNe spectra show stronger Ha emis- 
sion than in TDE2. Furthermore, the spectrum of TDE2 
was taken at least ~ 40 days after the peak - if the peak 
occurs between the CSS and CRTS observations - and 
otherwise the spectrum is at least 100 days post-peak, 
yet the blue continuum shape seen in spectrum of TDE2 
is observed for only tens of days afte r the p eak for the 
type Iln SN in sample of IKiewe et"all (|2010f) . 

Comparison of TDEs to exotic SNe — We examined above 
the SNe Iln in the small but unbiased IKiewe et al.l (|2010l ) 
sample. Here we consider specific SNe reported in the lit- 
erature; due to selection biases these tend to be particu- 
larly luminous and unusual. Only exceptionally blue and 
slowly decaying, relatively luminous SNe could mimic 
our TDEs. An exa mple of a slow ly decaying and lu- 
minous SN is 2006gy (|Quimbvll200(l . which' was a rather 
bright explosion (Mr = — 21.3±0.1) tha t is considered to 
be an interaction-dominated type II SN (|Agnoletto et al.l 
120091 ). Immediately after its peak, the decay of 2006gy, 
at 0.2 mag day -1 , is much faster than in either of our 
TDEs (see Fig. EQ]). At 270 days after its peak, the de- 
cay rate of 2006gy slows down to 0.004 mag day -1 , which 
is similar to known SN e of type II such as 1995G (Iln) 
(jPastorello et al.l 12002ft and 1999E (Ha, hybrid class) 
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(jRigon et al.l l2003h . This might lead one to speculate 
that our TDEs can be explained as emission from the 
late-time tail of the light curve of type II SNe. However, 
by the time their decay of flux becomes similar to our 
TDEs', the colors of SN 1995G, 1999E and 2006gy are 
red: B - V > 0.8. This is ve ry different from T DE1,2 
whose colors can be converted pester et al.ll2005l) to find 
B — V = —0.12, —0.02. Furthermore, the spectrum of 
2006gy at t > 270 days is characterized by large and 
broad emission lines, which is inconsistent with the flare 
spectrum of TDE2 (Fig. ITS]) . 

An example of a blue and UV brig ht SN is 2008cs 
(|Gezari et al.l l2009at iMiller et all l2009h . However the 
fast cooling of this SN implies a mean color of B—V = 0.6 
over the 60 day period after the peak, which is far more 
red than our TDEs averaged over a similar period. By 
day 91 the UV had faded 5 magnitudes fainter than the 
optical peak, in contrast to the ~ 2 magnitude decrement 
seen in the TDEs after ~ 800, 400 days. W e can also com- 
pare the TDEs to the transient SC P 06F6 (iBarbarv et~al] 
120090 . This has been suggested bv lQuimbv et al.l (|2011f ) 
to be a high-rcdshift example of a new class of blue 
transients that mark the deaths of the most massive 
stars. Using z = 1.19 for the redshift of SCP 06F6 
(|Quimbv et al.ll20Tll ). we extract a rest-frame u-band de- 
cay rate of 0.08 mag day" . This is an order of magni- 
tude faster than the decay rate of TDE1 and TDE2 in 
their rest-frame. Furthermore, the proposed new class 
of SNe should origin ate from faint, metal-poor galaxies 
(|Quimbv et al.l 120 111) , which are very different from the 
hosts of TDE1 and TDE2. 

Finally, we consider SN 2003ma ([Rest et al.ll2009l ). an 
extreme type Iln SN. The peak magnitude of the SN, 
Mr = —21.5, is similar to that observed in TDE2, 
and its decay rate during the first 300 days after the 
peak (~ 0.008 mag day -1 in the Ui?sM band, centered 
at 6250A) is also similar to TDE1,2. The SN was de- 
tected at d « 460 ± 85 pc projected distance from its 
host. Like the TDE2 spectrum, the spectrum of 2003ma 
shows no emission and no P-Cygni absorption up to the 
last detection at 1423 days after the peak of the SNe. 
However the intermediate width Ha emission is much 
stronger than what is observed for TDE2. Another dif- 
ference is that SN 2003ma is redder than TDE1,2: its 
color averaged over the first 50-100 days after the peak 
is B - R « 0.35 compared to B - R = -0.35, -0.18 
for TDE1,2. The SN also shows strong evolution of the 
foreground-extinction-corrected colors: V—I w at peak 
to V — I > 1 one year after the peak. There were no UV 
observations of SN2003ma but already by day 213 the 
B magnitude was 3.5 magnitudes below the peak in the 
I band, suggesting that by days 400-800 no UV would 
be seen. We not e that the host g alaxy of 2003ma is a 
starburst galaxy (|Rest et al.l I2009T) . as expected for the 
presumably rather massive (M > 10M Q ) progenitor of 
the SN; the star formation rate, derived from the nar- 
row line Ha emission of the host of SN 2003ma is one 
(five) orders of magnitude greater than the (upper limit 
on) the star formation rate from the TDE2 (TDE1) host 
spectrum. 

We conclude that TDE1,2 are not members of a known 
type of exotic SNe. However new types of SNe are dis- 
covered on a yearly basis so we have to consider the pos- 



sibility that we discovered a new exotic class of stellar 
explosion. In the following section we will use the dis- 
tance between the host and the flare to evaluate this 
possibility. 

4.2.4. Rejection of known SNe 

Geometric Rejection of stellar- distributed objects in nuclear 
sample — Although we rule out that TDE1,2 are ordinary 
SNe on the basis of the properties of the flares, it is of in- 
terest to know the chance probability for two flares whose 
progenitors follow the stellar distribution to be found as 
close to the nucleus as observed. We assume that the rate 
of ordinary SNe is proportional to the stellar light. This 
is justified by the good quality of the fit to the host-flare 
distance distribution (Fi g. \5j and the existi ng litera- 
ture on SN distributions. iFruchter et al.l (|2006f ) conclude 
that, while gamma-ray bursts are more concentrated on 
the brightest regions, the probability of a SN type la ex- 
ploding in a particular pixel is roughly proportional to 
the surface brightness of the galaxy at tha t pixel. This re- 
sult was confirmed bv lKellv et all ([2001 using also type 
II SNe; they conclude that both type la and type II SNe 
follow the galaxy light measured in the (/-band, with a 
clear exception being the rarer SNe Ic associated with 
long-duration gamma-ray bursts. (The observation by 
lAnderson fc Ja mes (20Q3) that SN type II do not trace 
star formation estimated from H Q + [Nil] emission is not 
relevant for this work, because our model relies on stel- 
lar lig ht, not star formation. Indeed lAnderson fc Jamesl 
(|2009| ) conclude that, except for a central deficit, type 
II SNe seem to follow the i?-band light, while SNIb/c 
appear more centrally concentrated with respect to the 
stellar light.) 

Taking the SN distribution to be given by the fit to 
the d distribution shown in Fig. [5] the expected con- 
tamination in the TDE sample is 0.9 SNe, in the mean. 
The probability that TDE1 and TDE2 arc ordinary SNe 
found by chance as close as observed to their hosts is 
the product of two factors. First, the Poisson probabil- 
ity of finding two or more SNe, in the TDE candidate 
(nuclear) sample, when 0.9 are expected; this probabil- 
ity is 23%. Furthermore and independently, if TDE1,2 
were ordinary SNe in the nuclear flare sample, the sepa- 
rations from the centers of their hosts would follow the 
distribution of stellar light within d < 0'.'2. The chance 
that two objects, drawn randomly from this distribution, 
have d values less than those of the TDEs, is w 7%. Thus 
the probability that the TDEs are actually ordinary SNe 
in the sample of flares passing our selection criteria, yet 
are so perfectly centered on their hosts as they are, is 
~ 1.6 %. This is the purely geometric suppression factor 
that must be applied to any hypothesized exotic type of 
flare whose progenitor follows the stellar light distribu- 
tion. 

Likelihood of SNe Iln — We have seen above in Section 
14.2.31 that individual cases can be found of SNe Iln which 
come close to matching particular properties of TDE1,2, 
but that when the ensemble of observations - luminos- 
ity, late time UV emission, spectrum, color evolution 
- are considered, TDE1,2 bear no resemblance to any 
SN observed to date. Nevertheless, the properties of ex- 
treme SNe depend strongly upon their environment, so 
one might think that the TDEs could simply be the lat- 
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est exotic SNe, whose properties differ from any SN seen 
earlier. In such a scenario, TDEI,2 would be on the tail 
of the distribution of SNe Iln, which are themselves a 
small fraction of all SNe. Only 17% of SNe in a flux lim- 
ited search such as ours are type II dLi et al.ll20loh . of 
which about 29% are type Iln and extreme examples are 
still more rare. The radial distribution of type II SNe 
is approximately the same as for all SNe, c.f., Fig. 8 in 
(|Leaman et al.l f2010h . so the 1.6% geometrical penalty 
for the SNe occurring so close to their galaxy centers ap- 
plies. Thus the probability of finding two type Iln SNe 
hosts and such central locations as TDE1,2 is < 0.08 %, 
and is still lower for more exotic types of SNe for which 
TDE1,2 might be first examples. 

We note that SN2006gy is a type Iln SN which would 
have passed our pipeline selection criteria had it been 
at the redshift of TDE1,2: its proximity to the center 
of its galaxy would have placed it in the nuclear sam- 
ple, being resolved only due to its lo w redshift and high - 
resolution imaging as at re 350 pc (jSmith et al.l 120071 ) . 
and its host would not have been excluded by the QSO 
locus. SN2006gy is readily recognized as a SN and dis- 
tinguished from TDE1,2 based on the color-evolution of 
the flare, as discussed in Section l4. 2.3. 01 above, but it is a 
concrete reminder that geometrical background rejection 
alone is not sufficient to eliminate SNe in a sufficiently 
large sample of TDE candidates, and in practice flare 
properties need to be considered. 

4.2.5. New SN type, only found near the centers of galaxies 

Since the properties of TDE1,2 are unlike any known 
SN, one could entertain the idea that there may be some 
special type of previously-undetected SNe which occur 
exclusively as close or closer to the center of galaxies 
than TDE1,2. In that case the above statistical argu- 
ments would not apply. To assess the viability of such an 
option, in this section we make the ansatz that TDE1,2 
are examples of a new class of centrally concentrated SNe 
and see if the properties of the required class are reason- 
able. The distance between TDE 1,2 and the center of 
its host is d = 0'.'058, 0'.'068 with a mean accuracy of 
re 0'.'06; this corresponds to 0.14, 0.26 kpc in projection 
(but consistent with zero), respectively. 

Nuclear star clusters (NCs) are a rathe r mysterious 
phenomenon, observed in all Hubblc-types ()B6kerll2010( ) 
and studied in con siderable detail in the Milky Way 
dGenzel et al.ll201C). Th e nuclear star clusters reported 
by iWalcher et all (|2005D in nine late-type spirals range 



from 8 x 10 5 - 6 x 1O 7 M with a median re 1O 6 5 M ; 
the mean age of the last star formation burst is 34 Myr 
with the yo ungest stellar populat ion having a mean mass 
of 10 5 - 5 M^ (|Walcher et al.ll2006h . Having a radius of re 
5 pc, NC's provide a concrete example of a possible sce- 
nario with concentrated sources. 

If every galaxy contains a nuclear star cluster, we can 
infer the minimum required rate of TDE-likc explosions 
as follows. The average stellar mass of the galaxies mon- 
itored in our search is re 10 10 Mq, so re 3 x 10 -4 of all 
stars could be in nuclear star clusters. The total number 
of SNe that would have been detected in the TDE search 
if we had not imposed a cut on the host-flare distance or 
additional quality cuts is re 1500 Thus if a new class of 

14 The total number of SNe that would be detected if the d < 1" 
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Fig. 21. — Comparison of the TDE2 flare spectrum (ho st sub- 
tracted) to the two type Iln SNe from Kicwe ct al. (2010) which 
are most similar to TDE1,2 in terms of cooling rate and color. The 
spectra are normalized to be equal at long wavelengths. The legend 
lists the time difference with respect to the peak of the light curve. 
For TDE2, the location of the peak in the light curve is unknown 
but we have a low er limit based on the first CSS detection of the 
flare (see Fig. 1101 and discussion in Section I3.3.3H . The bump at 
H a in the TDE2 spectrum can be fit by a Gaussian with FWHM 
~ 8 X 10' 5 kms . The blue continuum of the prc-peak spectrum 
of SN 2005cp is similar to that of TDE2, but this SN spectrum 
shows more narrow and much stronger Ha emission. The P-Cygni 
absorption profiles, a classic feature of SNe spectra, are clearly 
inconsistent with the TDE2 spectrum. 



"nuclear" SNe is to explain the existence of two SNe we 
require the average rate per star of TDE-like SNe in a 
nuclear star cluster to be a factor n larger than the rate 
of normal SNe per star, where 



k 3 x 10" 



2/150 



(5) 



giving k re 45. 

This factor, k = 45, is the enhancement factor by 
which stars in the nuclear star cluster must explode as 
a SN of the new type, compared to having a normal SN 
explosion, if the TDE candidates are to be explained as 
first cases of a new class of nuclear SNe. Much more chal- 
lenging theoretically, is to explain the absence of TDE- 
like SNe at larger radii, where none are observed. With 
re 3000 times as many stars in an average galaxy as in its 
nuclear star cluster, we would expect to have seen re 6000 
TDE-likc events at larger radii if the rate per star were 
the same, requiring a suppression in the rate-per-star by 
a factor of at least 2.3/6000 compared to the rate in the 
nuclear star cluster (2.3 being the 90% CL upper limit 
on the number if are seen). 

Although NC's contain on average younger, more mas- 
sive stars than in the ensemble of monitored galaxies, a 
mechanism which produces a factor ~ 60,000 contrast 
between the rates of this new type of SNe in stars in a 
nuclear star cluster, compared to their rate elsewhere in 
the galaxy, may be difficult to devise. We note that the 
unusual SN 2006gy, w hich occurred at a bout 350 pc from 
the center of its host (|Smith et al.ll2007j) , could be an ex- 
ample of a "galatic nucleus SN" , but a larger sample of 

requirement were not imposed can be estimated by integrating the 
fit to the d-distribution (Fig. to d — > oo and taking into account 
number of SN lost by requiring two detections after peak of the flare 
(in the g and r band only, to avoid a color bias) and restricting to 
host outside the QSO locus. 
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similar SNe is needed to asses this possibility. 

4.2.6. Summary of comparison of TDE1,2 to SNe 

We have compared TDE1 and TDE2 to "ordinary" 
SNe, to an unbiased sample of type Iln SNe, and to par- 
ticular exotic SNe. Ordinary SNe have entirely different 
properties than the TDEs in their color, color evolution 
and decay rates. Although some SNe Iln and exotic SNe 
resemble the TDEs in some respects, e.g., rate of change 
in magnitude and color, there are no examples of SNe 
whose properties are not significantly different than the 
TDEs in at least some important respects. Flares from 
known SNe types including Iln's follow the stellar distri- 
bution to a good approximation, with the result that the 
probability of two SNe Iln flares occurring so close to the 
nucleus as TDE1 and TDE2, is < 0.08 %. 

Thus if TDE1,2 are SNe, they must be the first exam- 
ples of a new exotic type of flare which occurs only very 
near the center of galaxies, to evade the 2 x 10 -2 geomet- 
rical penalty and explain the failure to detect many more 
events at larger radii. Since known SNe are essentially lo- 
cal phenomenon - depending only on the progenitor star 
and its immediate environment - and there are so many 
more stars in the bulk of the galaxy than in any physical 
system near the galactic center, a model along these lines 
requires that either the progenitors of TDE-like SNe are 
concentrated by a large factor near the nucleus of the 
galaxy or some process in the nucleus stimulates TDE- 
like explosions at a much higher rate than in ordinary 
SNe. We considered nuclear star clusters as a possible 
example, and found a minimum required enhancement 
factor of « 45. Equally or more challenging, the rate of 
these new SNe must be at least a factor w 3000 lower for 
stars outside the nuclear star cluster than for stars in it. 

Given the evidence above, we conclude that TDE1 and 
TDE2 - the two events which survived our pipeline cuts 
- are unlikely to be supernovae. 



5. DISCUSSION 

In this section we first compare the observed properties 
of TDE1,2 with properties of TDE candidates in the lit- 
erature, and then compare with theoretical predictions. 

5.1. Comparison to TDE candidates in the literature 

IGezari et al.l (|2009b[ ) presented two candidate TDEs 
discovered with GALEX that had simultaneous weekly 
optical difference imaging. The optical colors and light 
curve shape of TDE1,2 are similar to the g, r, and i and 
NUV light curves of those two GALEX TDE candidates. 
The black body temperatures of those GALEX candi- 
dates, measured by fitting to both UV and optical data, 
are ~ 5 x 10 4 K. This is to be compared to - 2 x 10 4 K for 
TDE1,2, obtained by fitting to optical data only. Thus 
the observed optical properties of our two TDEs are sim- 
ilar to each other and to the GALEX TDE candidates 
Given the small number of examples, the differences can 
be attributed either to variations within the TDE popu- 
lation or to the differences in the selection biases of the 
observations. 

Unfortunately, simultaneous optical imaging is not 
available for the X-ray selected TDE candidates, hence 
a comparison is possible only for derived properties. As 



discussed in IGezari et al.l (|2009b[) , the black body tem- 
peratures derived from the X-ray spectra are higher than 
those obtained for the UV and optical-selected TDE can- 
didates, which can be explained by the geometry of the 
emitting region (i.e., the high energy photons are pro- 
duced closer to the black hole). 

5.2. Comparison of flare properties to theory 

Predicting the properties of optical TDE flares is non- 
trivial, and the range of variation in observed flare prop- 
erties for a given black hole mass is expected to be 
large due to the variety of possible pericentric distances, 
masses and radii of the disrupted star, black hole spin, 
and viewing angles. Therefore, a much larger sample of 
TDEs will be needed - selected as done here based on 
host and geometrical criteria rather than flare proper- 
ties, with spectra and much more finely sampled light 
curves and simultaneous observations outside the optical 
band - before the full complexity of the tidal disruption 
phenomenon can be even superficially explored. 

Here we give a preliminary comparison of the most ba- 
sic properties of the flares - summarized in Table [2] - to 
the pr edictions of t wo theoretical models, lLoeb fc Ulmerl 
(fl997l LU97) and IStrubbe fc Quataerti (120091 . SQ09T 
LU97 is a simple model based on thermal emission at the 
Eddington luminosity, while SQ09 follows the evolution 
of the tidal debris and makes detailed numerical predic- 
tions for emission during and after the super-Eddington 
period, in three illustrative examples. 

The most massive black hole in the SQ09 examples has 
Mbr = 10 7 Mq - in the middle of the estimated black 
hole mass range for the TDE1 host galaxy, but at the 
very low end of the range for TDE2. For this example, 
SQ09 predict a peak value of vL v ss few x 10 41 erg/s in 
the g-band. This is 1 (2) orders of magnitude below the 
observed, i.e., post-peak maxima of TDE1 and TDE2, 
respectively. Furthermore, SQ09 predicts a temperature 
which is considerably larger than found with a black body 
fit to the optical SED. 

The discrepancy in the SQ09 modeling is particularly 
severe for TDE2 because i) according to SQ09 Fig. 4, 
the peak luminosity decreases for larger Mbh and the 
central mass estimate for the TDE2 black hole is higher 
than in the SQ09 example, ii) we know from the CSS 
prc-SDSS observation that the flare is at least w 40 and 
possibly greater than 90 days old when seen by SDSS, de- 
pending on whether the peak occurs after or before the 
CSS observation, and iii) the observed g-band luminosity 
is a factor « 100 higher than predicted at the peak for 
this example. (Interestingly, the peak optical luminosi- 
ties and blackbody tem peratures of t he TDE candidates 
discovered by GALEX (jCezari et ai1l2009bL Table 2) are 
also inconsistent with SQ09 predictions.) 

Table [3] reveals that the color and time evolution of 
TDE1 and TDE2 are quite similar to each other. The 
pipeline selection process did not place constraints on 
the color and, although too-short flares would not be 
accepted, the lifetimes of TDE1,2 arc far longer than re- 
quired to pass the selection criteria. Thus it is reasonable 
to think that these properties are typical for optically dis- 
covered TDE flares. Luminosities on the other hand are 
naturally biased toward the high end of the distribution. 
Since the volume of detectability of a flare increases as 
L 3 / 2 , we cannot exclude the existence of a population of 
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dimmer flares. 

This raises the question of how to account for the ob- 
served flare properties. The earlier and simpler model 
of LU97, based on thermal emission at the Eddington 
luminosity, correctly predicts the observed temperature 
(within uncertainties) and its slow evolution, but pre- 
dicts a much higher luminosity than observed. LU97 
also predicts that the luminosity decays much more 
slowly than observed. The LU97 model can be rec- 
onciled with observation if their hypothesized optically 
thick envelope subtends only a fraction of the An solid 
angle, and this fraction decreases with time. Or per- 
haps SQ09 arc on the right track but they adopted 
parameter choices which need modification, or bet- 
ter modeling of the r adiatively driven wind is needed 
(jLodato fc Rossil 120111 ). Or some important aspect of 
the process may be missed completely; for instance the 
presence of a weak, pre-existing accretion disk might sig- 
ni ficantly enhance the powe r of the flare, as proposed 
bv iFarrar fc Gruzinovi (|2009l) as an explanation for the 
correlation repor ted between ultra-high energy cosmic 
rays and AGNs (jThe Pierre Auger Collaboration! 120071 
120081 ) whose luminosities are t oo low to accelerate pro- 
tons to the obse rved energies (jFarrar fc Gruzinovi [200l 
iZaw et al.ll2009l) . 

Note added: After the properties of TDE1,2 were pre- 
sented in the prelimin a ry archive version of this w ork 
(jvan Velzen et alJl2010T ). iStrubbe fc Quataertl (|2011f ) ex- 
plored a larger parameter space in their model and found 
that TDE1 can be readily explained adopting adjusted 
parameter assumptions. They report that the luminosity 
of TDE2 can be explained in their framework as well, but 
agree with our observation that its color and slow decay 
are more suggestive of the state considered by LU97. 

The luminosity and temperature of a tidal disruption 
flare, and their time evolution, depend on the (unknown) 
time between disruption and first observation, and on 
unknown or poorly-known parameters of the black hole 
(-Mbh and spin) and the initial star and the orbit and 
viewing angle. This means that the range of flare types 
is enormous and one might despair of being able to 
test theoretical models. Remarkably, however, the ra- 
tio of the cooling rate to the bolometric luminosity de- 
cay, (d\nT/dt)/(d\nL /dt), is independent of all these 
unknown parameters and also independent of the time 
since disruption, in the SQ09 model during the super- 
Eddington phase. Combining equations 2, 11, and 13 of 
SQ09 we derive that 

(d\nT/dt)/(d\nL/dt) = -5/4 (6) 

for any TDE, at any time during the super-Eddington 
phase. This is certainly not the evolution observed 
for TDE1,2 in the g-band, for which we measure 
(dlnT/dt)/(dlnL q /dt) = ± 0.2 and 0.2 - 0.7 respec- 
tively. Whether this is further evidence of a problem with 
the SQ09 model requires detailed modeling since the evo- 
lution of the <?-band luminosity may not follow that of 
the bolometric luminosity, and a single black-body may 
not correctly describe the SED. We present the relation 
Eq. |6]here because of its power to test the picture of the 
wind-driven super-Eddington phase independently of the 
initial conditions of the disruption event. Simultaneous 
measurement of a larger portion of the SED, to allow 
Eq. © to be tested, would therefore be highly beneficial 



in future observational campaigns to explore the TDE 
phenomenon. 



6. FUTURE SURVEYS 

We estimate below the detection rate of TDEs which 
can be expected in current and future optical transient 
surveys, for the pipeline used here and similar observa- 
tional conditions and cadence as SDSS Strip e 82. This es- 
timate differs from earlier estimate s such as (jGezari et al.1 
I2009bt [Striibbe fc Quataertl 120091 ) because those ignore 
the cost in event-rate implied by cuts needed to insure 
clean and unambiguous detections, and our estimate here 
incorporates the observed properties of TDEs rather than 
relying on models. 

The total effective time spanned by the SDSS Stripe 
82 data, r b s , is just the total time between observations 
within a season; the > 9 month gap between seasons 
is not included in r b s - To account for the difference 
in sampling across Stripe 82 we calculate r b s in bins 
of width 3.6 degree along right ascension. The mean 
T bs is 1.03 yr with a mean cadence of 7.5 days. Thus 
two detected TDEs corresponds to a TDE detection rate 
of A^obs = 1.9 yr -1 . We can scale this detection rate to 
current and future optical surveys of similar cadence and 
selection criteria using the flux limit, Fy lm , and faction of 
the sky observed, / s ky: 

Aobsa/skyi^ 72 • (7) 

This yields A obs = 13, 14, 4180 yr" 1 for the PTF, Pan- 
STARRS Medium Deep Survey and LSST, respectively. 
In Table [5] we list the adopted values of Fy im , / s ky These 
detection rates are lower bounds on the actual number of 
TDEs which can be observed (if the observational quality 
is equal to that of SDSS and these TDEs have typical 
luminosities) , because the cadence can be optimized and 
the pipeline made more efficient, to maximize detections 
for any targeted light-curve type in a dedicated survey. 

Future optical surveys will be predominantly photo- 
metric and will generally not have a large fraction of 
hosts for which spectra have been obtained, as we have 
for SDSS Stripe 82. In fact, this need not prevent ob- 
taining a TDE candidate sample for followup with 0(1) 
false positives, if the angular resolution allows adequate 
rejection of non-nuclear flares. The first line of defense 
against SNe contamination is good resolution. The pu- 
rity of the nuclear sample is determined by the accuracy 
with which the flare-host separation can be measured, 
because SN background increases very rapidly as the res- 
olution is compromised. 

When the goal is rapid, intensive spectral and multi- 
wavelength follow-up rather than discovery in archival 
data, the appropriate selection strategy changes from the 
one used for here for TDE1,2. The first priority is to be 
confident that a flare passing the selection criteria has a 
very high chance of being a TDE and low chance of being 
uninteresting. Elements of such a strategy are: 
• QSOs can be suppressed in the target sample by ex- 
cluding galaxies within the QSO locus in Fig. [H] galaxies 
showing continuing irregular variability during monitor- 
ing can be excluded as presumptive AGNs without spec- 
troscopic follow-up, if spectroscopic resources are limited. 
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TABLE 5 

Detection rates of other optional surveys. 



Survey Name 


Cadence 




/sky 


N ohs 






(mag) 


(yr- 1 ) 


CSS [1] 


14 days 


19.5 


0.6 


5 


MLS [1] 


14 days 


21.5 


0.09 


12 


QUEST [2] 


hours to years 


20.5 


0.36 


12 


Palomar Transient Factory [3] 


5 days 


21.0 


0.2 


13 


Pan-STARRS Medium Deep Survey [4] 


4 days 


24.8 


0.0012 


15 


Pan-STARRS 3tt Survey [4] 


6 months 


23.5 


0.75 


1557 


Large Synoptic Survey Telescope [5] 


3 days 


24.5 


0.5 


4131 



Ref ere nces. — [1] (Drake et all d2009t ). [2l lHadiivska et al.l ll2TjTlh . [3l lLaw et all ||2009|) . [41 IGezari et~aT1 
( 2008). r5l llvezic et all 120081 



Note. — The survey plus reference used to obtain or estimate the cadence, flux limit (-Fn m ) and fraction 
of the sky covered (/ s ky) are listed. We scale the detection rate using, Eq. [7] and N a t, s = 1.9 yr for the 
analysis presented here. We have used 300 deg 2 as the angular area for S tripe 82. Since the cadence of the 
observations of Stripe 82 decreases towards the edges, ISesar et all 120071 ') have used 290 deg 2 for this arc;). 
However the total area of Stripe 82 that is imaged is 312 deg 2 ; we thus adopted 300 deg 2 as a reasonable 
value to obtain f s ^ y . 



• TDE1,2 both fall in a"TDE-locus" based on photo- 
metric properties of the flares alone (Fig. [T2"|) . allowing 
powerful rejection of AGN flares and SNe without spec- 
troscopic followup. Fig. [20] shows that flares from the 
hosts which are identified as QSOs based on their pho- 
tometric properties only, have properties similar to the 
flares from spectroscopically confirmed QSOs. Moreover 
requiring a flare to fall in a "TDE-locus" of very blue 
flares - u- g < -0.1 and g-r < -0.2 (Fig. [12J - re- 
duces the contamination of SNe in a TDE search by a 
factor of > 50. Thus photometry alone, without prior 
spectroscopy, is sufficient to reduce the variable-AGN 
contamination of the TDE candidate sample, and rejects 
almost all SNe. 

• The only SNe which are observed in the TDE-locus 
in flare color are of type Iln, but these occur only very 
rarely in E/S0 galaxies (|Li et al.l [2010h . Thus without 
introducing a bias in the selection of TDEs, targeting 
can be restricted to early- type galaxies. This reduces 
the contamination of SNe Iln, the most troublesome SN 
contaminant in a TDE search, based on color alone. The 
appropriate trade-off between stringency of rejection of 
SNIIn and loss of real TDEs will be determined by the 
specifics of the survey and follow-up resources. 

7. SUMMARY 

We have presented here the first two compelling, 
optically-discovered, stellar tidal disruption candidates. 
This work demonstrates the feasibility of discovering 
TDEs in optical synoptic surveys without imposing se- 
lection criteria depending on the properties of the flare. 

Our pipeline rejection of non-TDEs is based on ge- 
ometrical criteria to eliminate supernovae and on host 
properties to eliminate variable AGN flares, rather than 
being based on properties of the flares themselves, in or- 
der to minimize selection bias. The pipeline rejection 
gives an a priori probability of the flares being SNe or 
variable AGNs of less than 3%. Thanks to the very large 
sample of galaxies in SDSS Stripe 82, with a large num- 
ber of them having spectra, we have excellent data on the 
properties of host galaxies and their variability. This al- 
lows us to remove 90% of the QSO hosts by a color-color 



cut on the host galaxy. The most serious remaining back- 
ground comes from variable emission from hidden active 
nuclei. These cases are excluded by the variability of 
their hosts in the non-peak seasons, eliminating 3 of the 
5 TDE candidates which survived the SNe geometrical 
and QSO color-locus cuts. Based on the variability ob- 
served in QSOs and identified Seyferts, we estimate the 
probability that a QSO or Seyfert which satisfies our flare 
selection criteria in some season, shows as little variabil- 
ity in the other seasons as TDE1,2, to be < 2 x 10 -5 . 

Further SDSS observations, follow-up spectra and 
GALEX observations provide powerful a posteriori ev- 
idence that TDE1,2 are neither SNe or variable AGNs. 
The host spectra are consistent with the hosts not hav- 
ing active nuclei. Although no requirement was placed 
on the flare properties in the selection process, the TDE 
flares are very distinctive in comparison to SNe and flares 
in variable QSOs and AGNs. The lower limit on their in- 
crease in flux compared to a possible AGN contribution 
to the baseline flux is far greater in the TDEs than any 
QSO flare in the Stripe 82 data (Fig. [18]). The mean 
color and cooling rate, as well as the decay rate, of the 
TDE flares arc significantly different from any SN in our 
sample, as shown in Figs. [T2land[T9l The TDE flares are 
significantly bluer than QSOs and QSO flares, although 
like QSO flares their color evolves very slowly compared 
to SNe. In particular, GALEX recorded a level of UV 
emission from TDE2 w 800 days after the flare that is or- 
ders of magnitude greater than in any known SN. GALEX 
observations of TDE2 combined with color information, 
show that TDE2's flare is unlike every known SN flare in 
at least one respect. The closest resemblance is to type 
Iln SNe and a few exotic SNe. 

Serendipitously, a spectrum of TDE2 was taken during 
the flare, a few days after the first SDSS detection of the 
flaring state. We have recently taken a spectrum of the 
host galaxy, giving us by subtraction a spectrum of the 
flare itself. Wc compared of this spectrum to spectra 
from an unbiased sample of type Iln SNe, finding that 
these SNe spectra show stronger Ha emission. 

The properties of the flares argue against the possi- 
bility that TDE 1,2 are an unusual but known type of 
SN. The most TDEl,2-like SNe are type Iln's, but the 
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probability of finding two of those, and in such central 
locations as TDE1,2, is < 0.08 %, and is still lower for 
more exotic types of SNe. A final option, that these are 
the first examples of a new class of SNe occurring only 
at the centers of galaxies is shown to require a thousand- 
fold or greater enhancement in the rate of such events 
in stars near the nucleus compared to the rate in stars 
located elsewhere. 

With only two examples of probable TDEs to study, 
it is still very much early-days for testing models in the 
literature. Nevertheless, the events have already enabled 
refinem ent of the recent detailed modeling of the pro- 
cess bv lStrubbe fc Quataertl (|2009l SQ09): the observed 
luminosities of TDE1,2 are at least 1-2 orders of magni- 
tude higher than predicted there, and the temperatures 
determined by a black body fit to the optical SED are 
considerably lower than in the simulations, for the as- 
sumed parameter choices. However with adjusted pa- 
rameter choices the features of TD E1 can be explained 
and t he luminosity of TDE2 as well (|Strubbe fe Quataertl 
1 2 1 If) . We show that, independently of initial conditions 
and the time since disruption, the SQ09 model in its 
bright, super-Eddington phase predicts that the temper- 
ature of a TDE flare increases while its bolometric lu- 
minosity decreases, with a specific relationship between 
the rates which is independent of initial conditions or 
time since flare. This cannot be tested directly for these 
TDEs because most of the luminosity is predicted to be 
emitted above optical frequencies, but when such mea- 
surements become available they will provide a decisive 
test of the SQ09 mod el. The earlier and simpler model of 
lLoeb fc Ulmerl (|1997|) . based on thermal emission at the 
Eddington luminosity, can be reconciled with observation 
if their hypothesized optically thick envelope subtends 
only a fraction of the 4ir solid angle, and this fraction 
decreases with time. 

Our work demonstrates that candidate tidal disrup- 
tion events can be efficiently identified using photomet- 
ric imaging alone. We conclude based on our observed 
TDE rate and pipeline efficiency, that current and next- 
generation optical synoptic surveys should contain hun- 
dreds or thousands of TDEs. We have shown that a 
TDE candidate sample with 0(1) purity can be identi- 
fied using host selection, geometric resolution and color 
alone. With such a sample, the cost of excluding im- 
posters with a follow-up observation the next night or 
later that night on another instrument is sufficiently low, 
that a campaign to create a large sample of tidal disrup- 
tion events with high-frequency, multi- wavelength obser- 
vations is feasible. This will allow the tidal disruption 
phenomenon to be explored in full detail, opening an ex- 
citing new chapter in black hole astrophysics. 
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